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Résumé 


Résumé 


Les interactions de la méthoxyphényl-N-sulfamoyloxazolidinone avec les ions Cu (I) et Co 
(D) dans l'éthanol à 25 ° C ont été étudiées expérimentalement par spectrophotométrie UV-Vis. 
Le traitement des données avec un ajustement des moindres carrées non linéaires a permis la 
détermination des stæchiométries, des constantes de stabilité des complexes et des diagrammes 
de distribution des espèces dans les milieux réactionnels en fonction de la concentration en 


ligand. 


Les complexes à l'état solide ont été synthétisés et caractérisés par FT-IR, RMN 'H, techniques 
d'analyse thermogravimétrique (TGA) et calorimétrique à balayage différentiel (DSC) afin de 


prendre en compte la décomposition thermique du complexe de Cu (I). 


L'étude DFT a été effectuée pour obtenir des informations sur l'interaction L: Cu (ID) et / ou Co 
(I). D'abord, les propriétés structurelles prédites du ligand ont été comparées à celles 
expérimentales obtenues à partir de données cristallographiques. Les performances de chaque 
fonctionnelle d'échange B3LYP, B3LYPD3, B2PLYP et B2PLYP-D3 avec la même base 6-311 
++ G (d, p) ont été testées dans la prévision des propriétés structurelles. Les paramètres de 


réactivité chimique globaux et locaux ont également été estimés. 


Les complexes métalliques les plus stables (1: 2) ont été optimisés au niveau B3LYP dans un 
milieu à l'éthanol avec le modèle CPCM. Les spectres FT-IR calculés et les résultats TD-DEFT 
sont en accord avec les données expérimentales correspondantes. Les interactions ont également 
été analysées et caractérisées par orbitales naturelles de liaison (NBO), potentiel électrostatique 


moléculaire (MEP) et non linéaire optique (NLO). 


Enfin, un docking (amarrage) moléculaire a été réalisé pour étudier les activités biologiques 


relatives du ligand seul et des complexes métalliques. 
Mots clés: 


N-sulfamoyl oxazolidinone, Complexes Cu (I) et Co (I), DRX, DFT, Docking moléculaire. 
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Abstract 


Abstract 


Interactions of methoxyphenyl N-sulfamoyloxazolidinone with Cu(I) and Co(Iions in 
ethanol at 25 °C were studied experimentally using UV-Vis spectrophotometry. The data 
processing with a nonlinear least square fitting allowed the determination of stoichiometries, 
stability constants of the complexes and species distribution diagrams for each complex against 
ligand concentration. 
The solid-state complexes were synthesized and characterized by FT-IR, THNMR, 
thermogravimetric analysis techniques (TGA) and differential scanning calorimetry (DSC) to 
account for the thermal decomposition of Cu(Il) complex. 
DFT study was performed to obtain insights on L: Cu(l) and / or Co) interaction. First, the 
predicted structural properties of the ligand were compared with the experimental ones obtained 
from crystallographic data. DFT calculations were performed at B3LYP, B3LYPD3,B2PLYP 
and B2PLYP-D3 levels with the same 6-311++G(d,p) basis set and consequently the 
performance of each exchange functional was tested in the structural properties prediction. Also, 
global and local chemical reactivity parameters were computed. 
The most stable metal complexes (1:2) were optimized at B3LYP level in ethanol medium with 
the CPCM model. Computed FT-IR spectra and TD-DFT results agree with the corresponding 
experimental data. The interactions were also analyzed and characterized by NBO, MEP and 
NLO. 

Finally, a molecular docking was performed to investigate the relative biological activities of 


the ligand alone and of the metal complexes. 
Key words: 


N-sulfamoyloxazolidinone, Cu(ID) and Co(Il) complexes, XRD, DET, Molecular docking. 
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La chimie est seule, peut-être, à avoir des frontières avec l’ensemble des autres disciplines. 
Cela lui donne la possibilité d’intervenir de façon pertinente et originale sur pratiquement tous 
les grands défis de l’humanité du XxI° siècle (alimentation, énergie, santé, environnement, etc.). 
En effet, si aujourd’hui les disciplines sont devenues essentiellement interdépendantes, la chimie 
est probablement celle qui a le plus souvent et le plus profondément multiplié ses incursions chez 
les autres — sciences de la vie et de la santé, physique et matériaux, sciences pour l’ingénieur, 
sciences de la terre et de l’environnement, qui toutes la sollicitent en permanence. Le 
développement de ces interfaces constitue sans nul doute l’un des enjeux les plus importants de 
la science contemporaine. Malheureusement, 1l n’est pas certain que la recherche scientifique de 
nos Jours soit à la hauteur de cette ambition tant les barrières culturelles et organisationnelles à 
la multidisciplinarité sont fortes. 

Sollicitée, la chimie l’est parce qu’elle est avant tout la science de la création des molécules et 
des matériaux dont dépend tout le reste. C’est donc elle qui façonne le monde concret dans 
lequel nous vivons (médicaments, cosmétiques, polymères, plastiques, verres, pour ne citer que 
quelques composés chimiques de notre univers). À l’évidence, la chimie du vivant, à la fois 
chimie bio-organique et chimie bioinorganique, jouera un rôle majeur dans ce nouveau contexte. 
Enfin la chimie biologique continuera à trouver des solutions originales dans le domaine de la 
santé. La biologie et ses applications en médecine ont constitué la grande aventure scientifique 
de la fin du XXe siècle. Le résultat est extraordinaire, comme le montre l’accroissement continu 
de l’espérance de vie, en particulier dans les pays développés. La chimie a joué dans ce progrès 
un rôle majeur. Il ne fait aucun doute que la recherche de nouveaux médicaments restera une 
activité centrale de la chimie du vivant dans les domaines de l’infectiologie, des grandes 
maladies neuro-dégénératives ou cardio-vasculaires et du cancer. Malgré les progrès des 
technologies de thérapie génique et cellulaire qui ouvrent des perspectives passionnantes, le 
médicament restera longtemps encore très majoritairement moléculaire — molécule de synthèse 
ou naturelle, petite molécule ou macromolécule biologique — et les progrès de la médecine et de 


la santé resteront largement tributaires de la chimie. 
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L'une des familles de molécules biologiquement actives est celle des sulfamides et de leurs 
dérivés. Particulièrement, les sulfonamides et leurs dérivés sont largement utilisés en médecines 
humaine et vétérinaire [1]. Ils ont été employés comme antimicrobiens [2], antifongiques [3], 
antipaludiques [4], activité antitumorale [5,6] et les inhibiteurs de l'anhydrase carbonique (sous 
forme de diurétique ou réactifs hypoglycémiques) [7-91]. Ils sont également utilisés comme 
agents pharmaceutiques pour le traitement de différentes maladies telles que les infections [10] la 
maladie d'Alzheimer [11] et le VIH [12]. 

Par ailleurs ces dernières années, en plus de l'importance des molécules libres, l’attention 
s’est focalisée sur l'implication de leurs complexes métalliques et organométalliques [19]. En 
particulier, les sulfonamides d'argent se sont révélés comme agents antimicrobiens topiques 
efficaces, notamment l'Ag-sulfadiazine utilisés dans le traitement des brûlures [20]. En plus, 
plusieurs séries de composés homoleptiques et hétéroleptiques de Cu (ID) ont été étudiées, pour 
montrer l’importance des sites de coordination pour l'activité antibactérienne [21,22]. Selon des 
travaux menés sur les complexes de sulfonamide Zn (I), Cu (D), Ni (D), Ce (D), Bi (D), Cd 
(D), Hg (D, SmD) et UO:() il a été trouvé que la polyvalence des sulfamides en tant que 
ligands et l’importance de leurs complexes dans la chimie de coordination et la chimie 
médicinale [23-28]. 

Dans le cadre de cette thématique, le présent travail se propose d'étudier les complexes de 
deux ligands : méthoxyphényl-N-sulfamoyloxazolidinone (L!) et N- sulfonyl imine (L?) avec 
les ions Cu (ID) et Co (ID. 

En solution d’éthanol et à 25 °C, l’interaction entre chaque ligand avec chaque ion a été 
étudiée par spectrophotométrie UV-visible. Le traitement des données par ajustement des 
moindres carrées non linéaires a permis la détermination des stæchiométries, des constantes de 
stabilité des complexes et des diagrammes de distribution des espèces pour chaque complexe en 
fonction de la concentration en ligand. 

Les complexes solides ont été synthétisés par FTIR, RMN 'H, techniques d'analyse 
thermogravimétrique (TGA) et calorimétrique à balayage différentiel (DSC) pour le suivi 
cinétique de décomposition thermique du complexe de Cu (I). 

L’aspect expérimental de ce sujet a été suivi par une étude complémentaire de modélisation 
moléculaire au niveau DFT. 

Le ligand Li a été optimisé aux niveaux B3LYP, B3LYP-D3, B2PLYP et B2PLYP-D3 sous la 
même base de 6-311++G(d,p). Les caractéristiques structurales des configurations obtenues 
(distances et angles) ont été comparées à celles obtenues de la DRX menée sur le cristal du 
ligand. Les résultats ont été discutés sur la base de la nature des fonctionnelle d’échange et de 


l’implication de la dispersion. Afin de prévoir les sites d’interaction avec les espèces 
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électrophiles ou nucléophiles, une prédiction des paramètres locaux de réactivité chimique a été 
menée. 

Les complexes les plus stables de stœchiométrie L'ion — 1:2, ont été optimisés au niveau 
B3LYP en milieu d’éthanol en adoptant le modèle CPCM. Les spectres FTIR et UV-visible des 
complexes prédits moyennant la DEFT et la TD-DET ont été comparés à ceux obtenus 
expérimentalement. Les interactions ont été également caractérisées et analysées par NBO, MEP 
et NLO. 

Finalement, l’étude d’amarrage (docking) a permis de déterminer l’orientation de la molécule 
du complexe étudié liée à protéines ciblés afin de calculer leur affinité et niveau d’activité. 

Le mémoire est présenté en quatre chapitres : 

Le chapitre I est constitué de généralités sur les complexes de coordination. Après 
présentation de la classification des ligands et des caractéristiques des métaux de transition dans 
les complexes de coordination, les diverses applications des complexes métalliques dans divers 
domaines notamment en biologie ont été exposées. 

Dans le chapitre IL, nous avons fait un rappel bibliographique sur les sulfonamides, les N- 
sulfamoyloxazolidinone et les N-sulfonyl imine utilisés dans ce travail, nous avons donné le 
mode de coordination de ces ligands dans ces complexes. Nous avons aussi, présenté une mise 
au point sur l’activité biologique essentielle des ligands et de leurs complexes. 

Dans le chapitre IIT, nous avons exposé les méthodes théoriques à savoir la modélisation 
moléculaire et le docking moléculaire utilisées pour étudier les complexes formés. Aussi nous 
avons présenté le principe général du docking moléculaire pour interpréter les interactions 
intermoléculaires ligand- récepteur protéiques. 


Le dernier chapitre est consacré aux résultats et discussion. 
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Les complexes métalliques et leurs activités biologiques 


La chimie de coordination est une discipline qui associe, le plus souvent, la chimie 
inorganique et la chimie organique. Plus précisément dans ce système, des molécules organiques 
appelées ligands sont liées à un atome ou ion inorganique central. Cette branche de la chimie a 
connu un développement important, non seulement dans le domaine de la chimie structurale 
(applications analytiques), mais également dans le domaine de la chimie pharmaceutique en 
raison des propriétés biologiques ou thérapeutiques d’un certain nombre de complexes. 

Grâce à leurs larges applications dans divers domaines d’activité chimique, biologique et 
environnementale, une attention particulière a été prêtée aux études des complexes de métaux de 
transition ces dernières années. 

Dans ce chapitre, nous avons décrit plusieurs travaux concernant les complexes de 
coordination, les types de complexes, la classification des ligands et l’application de complexes 


de coordination en biologie. 
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L.1.Historique de la chimie de coordination 


Par définition, la chimie de coordination est l’étude de composés chimiques résultant de 
l’interaction de fonctionnalités organiques ou inorganiques avec des centres métalliques. De 
telles fonctionnalités entourant l'atome de métal central sont appelées ligands à lier (à partir du 
verbe latin ligare). La synthèse et la caractérisation de tels composés de coordination sont 
essentielles à la compréhension des propriétés chimiques et physiques observées. 

A partir de 1893 et pendant les vingt années suivantes, Werner [1b, 2, 3] développa les idées 
modernes de la chimie inorganique. 11 a découvert que l'ion métallique (Co**) était entouré dans 
l'hexa amine cobalt de six ligands (6NH:), occupant une géométrie octaédrique. 

Alfred Werner était le pionnier de la chimie de la coordination (prix Nobel de chimie en 
1913[1b]). Ce savant a, notamment, étudié plusieurs complexes du cobalt avec des ligands amine 
et chloro, qui contenaient tous trois ions chlorure par ion cobalt, mais présentaient des couleurs 
et des propriétés physiques différentes (Figure 1.1). Ses observations ont conduit à la formulation 
que les atomes (ou molécules) peuvent être liés à un atome central dans certains motifs 
géométriques de la sphère de coordination primaire, tandis que les atomes de la sphère de 
coordination externe (ou secondaire) équilibrent la charge du complexe. Il a noté que le nombre 
d'atomes dans la sphère de coordination primaire, qu'il a appelé nombre de coordination, dicte 
souvent la géométrie autour de l'atome central. Par exemple, un nombre de coordination de six, 
noté comme étant particulièrement courant, donne généralement un complexe octaédrique, dans 
lequel les ligands de coordination occupent les sommets d’un octaèdre autour de l’atome central. 


Cela devint la base de sa théorie bien connue des composés de coordination. 


NH. 3 Cr ci 2Cr 
ax gx 
H3N NH; HN NH; 

3 H3 
| CH; 3 Cr ci Cr 

eg RQ 
H3N NH; H3N NH; 
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Figure 1.1 : Structures originales et révisées des complexes de cobalt étudiées par Werner 
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Aux composés de coordination, une grande importance a été accordée dans divers domaines 
de la science. Par exemple, le cis-platine, un complexe de platine intensément étudié depuis sa 
première synthèse rapportée en 1845 par Michele Peyrone (Figure 1.1). Ce complexe jouerait un 
rôle central dans le travail de Werner qui a correctement proposé la géométrie carrée et planaire. 
Il a été, également, en mesure de distinguer les isomères ci1s et trans, qui sont maintenant connus 
sous le nom de cis-platine et de trans-platine, respectivement [1,4, 5]. Actuellement, l'un des 
agents anticancéreux le plus efficace de la clinique, l'activité biologique du cis-platine, est dû à 
l'inhibition de la division cellulaire [6,7]. Le succès de la cis-platine a inspiré de nombreux 
chercheurs à mener des travaux visant à étudier le potentiel des complexes métalliques chimio- 


thérapeutiques [8-11]. 


I.2. Généralités 


Un complexe est un édifice polyatomique, constitué d’un cation métallique central (moins 
fréquemment d'un atome métallique), entouré d’ions ou de molécules associés à l’atome central 
par des liaisons chimiques. Le terme composé de coordination est, aussi, utilisé pour caractériser 
les complexes dont certains complexes sont constitués de plusieurs atomes centraux, appelés 


complexes polynucléaires. L’ensemble des complexes ont une formule générale comme suit : 
[MnLiXx] 


M: atome central (métal). 

L : ligands pairs (apporte une ou plusieurs paires d’électrons au métal). 

X : ligands radicalaires (apporte une ou un nombre impair d’électrons au métal). 

z: charge du complexe (si elle n'est pas nulle). 

I: nombre de ligands de type L coordonnés. 

x: nombre de ligands de type X coordonnés. 

La sphère de coordination est constituée par les ligands tandis que la sphère de complexation 

contient les ions complémentaires du complexe attirés par la force électrostatique. Autour de la 
sphère (autour de quoi), le solvant forme la sphère de solvatation qui stabilise et isole le 


complexe (Figure 1.2). 
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Figure 1.2 : Sphère de coordination 


Dans la formule d’un complexe, on indique l’ion métallique en premier, suivi des ligands 


chargés négativement puis les neutres et enfin ceux chargés positivement. 


L.3.Types de complexes 


La classification des complexes se base sur le nombre d’ions (ou d’atomes) centraux qu’ils 
comportent. Les complexes, dont les formules sont données ci-dessus, sont organisés autour d’un 
seul ion central. Ce sont des complexes monométalliques (on dit aussi mononucléaires). Si 
l’entité complexe comporte deux ou plusieurs ions métalliques, on ladésigne par les termes 


bimétalliques (binucléaire), trimétallique (trinucléaire), polymétallique (polynucléaire). 
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Figure 1.3 : Exemple d’un complexe polyatomique 


Dans les exemples ci-dessus, les ions métalliques sont éloignés les uns des autres. Lorsque la 
distance est faible, il peut se former des liaisons métal-métal et le complexe résultant (Figure 1.4) 


est appelé agrégat (cluster en anglais). 
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Figure 1.4 : Cluster bimétallique [Re2Cls]? 


L.4. Les ligands 


1.4.1. Généralités 


Un agent complexant (ligand) peut être défini comme une molécule ou un ion qui se 
coordonne à un centre organisateur qui est en général un atome ou un ion central. Ce ligand 
possède un ou plusieurs sites donneurs qui ne sont pas nécessairement tous engagés dans des 
liaisons de coordination. Ces ligands sont engagés par un ou deux électrons dans l’échange avec 
l’atome central. Il existe des ligands monodentés et des ligands polydentés selon le nombre des 
atomes liés avec l’atome central. Un halogénure est monodenté et l’éthylène diamine est bidenté. 


Un même ligand peut être coordonné de façon variable avec un atome central. 


1.4.2. Classification des ligands 


On peut classer les ligands selon le nombre d’électrons fournis sur le métal ou selon le 
nombre des liaisons formées avec l’ion métallique. Selon le nombre de liaisons que forme un 
ligand avec le métal, nous distinguons les ligands unidentés ayant une liaison avec le centre 
métallique (exemple: l’ammoniac NH) et les ligands polydentés ayant plusieurs liaisons avec le 


centre métallique. 


L.4.2.1. Les ligands chélatants 


Un ligand chélatant est un coordinat dont la structure est constituée de plusieurs atomes 
électrodonneurs, arrangés de manière à occuper plusieurs sites de coordination d'un seul ion 
métallique. L’exemple de la Figure 1.5 illustre le mode de coordination du ligand 1,2- 
diaminoéthane (= éthylènediamine, abréviation en), où les deux atomes d'azote sont liés au 
même métal M, formant un cycle stable à cinq chaînons. 

Du fait de la formation des cycles de chélation à cinq ou six chaînons, la présence de ligands 
chélatants dans les complexes métalliques engendre une plus grande stabilité comparée à leurs 


analogues monodentés [12]. 


10 


Chapitre I Les complexes métalliques et leurs activités biologiques 


Figure 1.5 : Complexe dérivé de l’éthylediamine, ligand chélatant 


L'exaltation de la stabilité d'un complexe par formation d'un ou de plusieurs cycles de 
chélation porte le nom d'« effet chélate », résultant essentiellement d'un effet entropique. Dans le 
cas du cation complexe [Ni(NB3) (0 )4 + (log B = 4,9), qui est nettement moins stable que le 
cation chélate [Ni(H>2O)4en]* (log B = 7,3), la formation du cation complexe provoque la 
disparition de deux espèces indépendantes. Alors que celle du cation chélate ne provoque la 
disparition que d'une seule espèce indépendante, induisant un accroissement relatif d'entropie 
pour la réaction de chélation, et donc un accroissement de stabilité [12]. Cette différence 
d'entropie, due aux différents nombres de particules impliquées et/ou aux constantes de vitesse 
plus élevées pour les cycles, serait due à l'effet de la grande concentration locale des 
intermédiaires monodentés [12]. 

Par ailleurs, une telle exaltation de la stabilité des complexes dérivés de ligands chélatants 
pourrait s'expliquer par des arguments microscopiques. En effet, un ligand chélatant agit comme 
une pince : si une extrémité se décoordine, l’autre extrémité la maintient proche du métal, la re- 
coordination est alors favorisée. De ce fait, le complexe issu d'un ligand chélatant polydenté, 


donc, tendance à moins se dissocier par rapport à son homologue. 


L.4.2.2. Les ligands pontants 


Un ligand pontant fait un pont entre plusieurs cations. Dans le di-u-hydroxobis(tétraaquofer 
()), aussi appelé octoaquo-di-u-hydroxo-difer(ID), (Figure 1.6),deux ligands hydroxyles font 
le pont entre deux cations métalliques. Au niveau de la nomenclature, la convention adoptée veut 
qu’un tel ligand soit précédé de la lettre «  ».Les ligands ponts les plus courants sont : OH, S?, 


NH2. 
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4+ 


Figure 1.6 : Complexe avec un ligand ponté (OH) 


L.4.2.3. Les ligands ambidentés : 


Un ligand ambidenté a deux atomes différents qui peuvent fonctionner comme ligand, c’est à 
dire que le ligand a deux atomes qui peuvent se lier avec le centre métallique mais leur 
arrangement stérique ne leur permet pas de former un anneau chélatant avec le métal. Ces 
ligands ambidentés peuvent donner lieu à des isomères de structure, citons comme exemple les 
ligands suivants : CN, CO, SCN, (CH3)SO (diméthyl sulfoxide=DMSO), HCON 
(CH3)2(diméthyl formamide= DMF). Ces ligands prennent un nom différent suivant l'atome par 
lequel ils sont liés. 

Aussi, 1l y a des ligands qui peuvent se lier au métal par un atome d’oxygène ou par l’atome 
d’azote. Ce ligand se nomme « nitro » s’il est lié par l’atome d’azote et & nitrito » s’il est lié par 


un atome d’oxygène (Figure 1.7). 


L #n-D- 
L. 46% 
L— di ” nitro 
| 
L 9 
L + 
| . _O| 
L— M— O—N7 nitrito 
| 
L 


Figure 1.7 : Complexe avec ligand ambidenté 


En effet, l’utilisation de ligands bidentés, possédant des hétéroatomes (N, O, S, P) combinés à 
un composé organique par exemple, a pour but d’améliorer la stabilité ainsi que les performances 
catalytiques des complexes [13]. Cela permet aussi de varier les effets de position des 
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groupements coordinnants sur les complexes. Les ligands bidentés conduisent généralement à 
une coordination en cis, tandis que celle de deux groupements monodentés sera essentiellement 
régie par leur encombrement respectif d'où leur effet trans. L'utilisation de ligands bi- ou 
tridentés permet également d’accroître la stabilité des complexes formés. En effet, la 
coordination d’un premier fragment du ligand faciliterait la coordination du second, par effet de 
proximité. Réciproquement, la décoordination du ligand serait plus difficile que dans le cas de 
ligands monodentés du fait de cet effet chélate [14]. 

Différents types de ligands chélatants polydentés ont été utilisés lors de la synthèse des 
complexes organométalliques et de leurs applications, aussi bien en catalyse, qu'en biologie, en 
magnétisme, etc. [15]. 

Ainsi, de nombreux systèmes organométalliques sont décrits et qui renferment, outre des 
ligands x donneurs, des substrats organiques chélatants à différents hétéroatomes, à l'instar des 
ligands de type N,N-donneurs (éthylène diamine, bipyridine, .)[16], des carbènes N- 
hétérocycliques [17], des ligands N, O, S ou P-donneurs bi- ou multidentés (bases de Schiff de 
type aniline salicylidènes et dérivés) [18], etc. Lorsqu'ils sont astucieusement combinés avec des 
ligands usuels, ils confèrent des propriétés catalytiques ciblées pour les systèmes du métal de 
transition résultant [19]. Ainsi, l'introduction de ligands chélatants de type base de Schiff par 
exemple (et notamment des aniline salicylidènes et dérivés) dans le complexe organométallique, 
devrait conférer des propriétés inédites à ce dernier et élargir davantage l'éventail de ses champs 
d'application [20]. 

Yudhvir et al. [21] ont aussi synthétisé et caractérisé un ligand base de Schiff tetradentate le 
bis (2, 5-dihydroxy acétophénone) éthylène diamine et ses complexes de Cu (Il) monomérique et 
dimérique (Figure 1.8). L'analyse par la spectroscopie UV-Visible montre une bande 
d'absorption du complexe entre 500 et 600 nm, due à la transition des électrons d'orbitales d-d. 
De plus, ils ont confirmé la coordination de Cu (IT) avec les sites d'oxygène (O) et d'azote (N) du 


ligand par la spectroscopie infrarouge (LR). 


ET Pa 
Cu 
X 5 LA 

HO o P OH HO (4 re) OH 

Cu A 
/I\ AN 

—=N — —N N— 

X X-Anion 


Figure 1.8 : Structures mono et dimérique des complexes synthétisés 
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L.4.2.4. Ligands macropolycycliques et cryptates 


Les ligands macropolycycliques [22] suscitent un grand intérêt car ils sont capablesde 
séquestrer différents substrats métalliques ou organiques pour former des complexesd'inclusion 
plus stables que les complexes macrocycliques. La majorité des composéspolycycliques décrits 
dans la littérature sont des macrobicycles ou des macrotricyclescomportant des atomes d'azote et 
d'oxygène ou de soufre. 

Ces cryptants ont unecavité tridimensionnelle permettant de fixer un ou plusieurs ions 
métalliques. Aïnsi, descomplexes macrobicycliques et macrotricycliques, mononucléaires et 
dinucléaires, ont étéobtenus avec les cations alcalins et les alcalino-terreux [23]. La stabilité de 
ces cryptatesest, généralement, élevée tandis que la sélectivité des récepteurs peut être ajustée 
eneffectuant des modifications structurales appropriées. Néanmoins, ces ligands présentent 
unefaible affinité vis-à-vis des métaux de transition par comparaison aux peraza-macrocycles 
[241]. 

De ce fait, la synthèse de ligands macropolycycliques polyazotés capables de séquestrer un ou 
plusieurs métaux de transition présente un intérêt majeur. En effet, de tels systèmes sont 
susceptibles d'être utilisés dans divers domaines tels que la catalyse, le transport ou la 
modélisation de sites actifs de métalloprotéines [25]. Dans cette optique, de nombreux ligands 
macrobicycliques et bismacrocycliques en série triaza- ou tétraaza-macrocyclique ontété 


synthétisés et complexés avec divers métaux de transition [25]. 


Figure 1.9 : Structure d’un macrotricycliques cylindriques 
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Figure 1.10 : Structure d’un macrotricycle sphérique 


L.4.2.5. Ligands triamines contenant des cycles azotés 


Les composés azotés hétérocycliques représentent une classe importante de ligands dans la 
chimie de coordination des métaux de transition (Figure 1.11). En particulier les ligands 
appartenant à la famille des azines (cycles azotés à six chaînons) tels que la 2,2°-bipyridine 
(bpy), la 1,10-phenanthroline (phen), la 2,2°,2”°-terpyridine (terpy) et leurs dérivés ont suscité 
depuis une cinquantaine d’années un immense intérêt aussi bien en chimie analytique qu’en 


chimie de synthèse organométallique. 


pyridine 1,10-phenanthroline 2,2°-bipyridine 2,2°,2°°-terpyridine 


Figure 1.11 : Exemple des composés azotés hétérocycliques 

Un tel engouement pour ces composés hétérocycliques polyazotés a tout naturellement orienté 
la chimie de coordination vers la synthèse et l’utilisation de ligands polyazotés contenant un ou 
plusieurs hétérocycles insaturés à cinq chaînons (azoles) (Figure 1.12). La nature des fragments 
hétérocycliques azotés coordonnés aux métaux de transition a une influence significative sur les 
propriétés chimiques, catalytiques et même physiques des complexes résultants. Cette influence 
n’est pas seulement liée au caractère s-donneur de l’atome d’azote présent dans les différents 
substituants azotés, mais dépend pour beaucoup du caractère p-accepteur/donneur du ligand. La 
grande majorité des composés de type azine sont d’excellents p-accepteurs pour les orbitales ‘d’ 
des métaux de transition, alors que les composés azoles sont plutôt de type donneur. Ces deux 
caractères p-accepteur/donneur donnent lieu à des transferts de charge très importants entre le 
métal et le ligand et influencent les propriétés d’oxydoréduction des complexes qui contiennent 
de tels ligands. 
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Un Ej EC 


pyrrole oxazole pyrazole imidazole 


Figure 1.12 : Exemple des composés hétérocycliques azotés 


Ainsi, vers la fin des années 80, l’étude de la chimie de coordination de ligands résultant de la 
combinaison de fragments azine/azole à été très importante, et de nombreux ligands bidentés N, 
N’-chélatants contenant des pentacycliques azotés ont été synthétisés et utilisés pour la chimie de 
coordination. Cependant, la synthèse de composés triamines tridentés constitués d’un ou de 
plusieurs substituants pentacycliques azotés a été moins développée et est encore maintenant 
restreinte à quelque composés-type. 

Les ligands triamines tridentés peuvent être classés en deux ensembles distincts, dont l’un 
englobe des composés de symétrie C2, tandis que l’autre contient des molécules tripodes 


triamines de symétrie C3. 


L5. Caractéristiques des métaux de transition dans les complexes de coordination 


Les complexes métalliques sont souvent classés par les ligands de coordination, qui peuvent 
contenir une variété d'atomes de coordination et de modes de liaison. L'état d'oxydation, qui 
indique le degré d’oxydation du métal, est au cœur de la réactivité du complexe métallique. Les 
deux électrons de la liaison métal-ligand étant attribués au ligand, de sorte que l'état d'oxydation 
correspond à la charge restant sur le centre métallique si tous les ligands étaient retirés. Les états 
d'oxydation formels sont souvent utiles pour classer les composés, ils ne sont pas des quantités 
mesurables car la plupart des liaisons comportent un certain degré de covalence, ou le partage 
d'électrons [26]. À cette fin, Jorgensen et Wieghardt [27] ont défendu l'idée d'une spectroscopie 
(ou l'état physique) d'oxydation, en particulier aux centres métalliques dans les complexes de 
coordination, qui sont mesurables et référencés à l'aide de données spectroscopiques et 
cristallographiques. Ce phénomène fait l’objet de nombreuses publications et il est connu pour 
des complexes présentant différents centres métalliques (3d* — 3d?) tels que les ions Fe(IID, 
Co(), Co(), Mn(), Mn(D) et Cr(I), l’ion Fe(Il) a fait l’objet de nombreuses études [28]. 

Les métaux de transition sont, par définition, des éléments à couches d incomplètes et leur 
dernière couche p (de valence) est vide. Pour donner lieu à l'existence de composés stables, ces 
éléments devront, plus ou moins, compléter ces sous-couches de valence par des électrons 


donnés ou partagés par les coordinats (ligands). Ces électrons apportés par les ligands leur 
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permettent ainsi d'atteindre, ou au moins d'avoisiner, la structure électronique du gaz rare suivant 
l'élément de transition dans la même période. 

Il est clair que l'environnement stéréochimique d'un ion métallique est un facteur important 
dans la détermination de son spectre. Un travail considérable a été accompli pour tenter de 
formuler des théories capables de décrire les liaisons dans les composés de coordination et de 
prédire leurs propriétés. Le premier succès, dans cette perspective, a été la théorie de la liaison de 
valence établie par Pauling et d'autres auteurs dans les années 30. Dans les années 50 et 60, la 
théorie du champ cristallin (TCC) et ses modifications, que l'on connait sous le nom dethéorie du 
champ des ligands(TCL) outhéorie du champ des coordinats, est devenueprédominante et a puisé 
son formalisme d'une autrethéorie qui est lathéorie des orbitalesmoléculaires. 

Bien que les théories de la liaison de valence etdu champ cristallin aient été 
largementdépassées comme modèles de travail, elles continuent à apporter leur contribution aux 
étudesactuelles sur les composés de coordination. Il faut donc savoir évaluer les modèles anciens 
pourbien comprendre les modèles modernes. 

Les quatre principales caractéristiques du métal de transition dans un complexe 
organométallique, permettant de bien définir formellement ce dernier sont : le Nombre 
d'Electrons de Valence (NEV), le Nombre d'Electrons Non Liants (NENL) et le Degré 
d'Oxydation (DO). 


L.5.1. Règle des dix-huit électrons 


En réalité, cette règle n’est autre que la règle de l’octet appliquée aux complexes de métaux 
de transition qui signifie que chaque atome souhaïitera mettre en commun ses électrons de 
valence pour atteindre la configuration du gaz rare qui le suit dans la classification périodique. 
Pour les éléments du bloc p, il suffit d’avoir 8 électrons pour satisfaire cette règle, la couche de 
valence étant complète pour ns’np‘.Alors que, pour les métaux de transition qui se caractérisent 
par la présence d’une sous couche d incomplète, cette règle classique valable en chimie 
organique est ici remplacée par la règle dite des 18 électrons car la sous couche d n’est complète 
que pour (n-1) d'ns'np°. Ainsi, en comptant les électrons «autour » du métal, on saura s’il peut 
ou non augmenter ou diminuer sa coordinence. La règle des 18 électrons n’est pas absolue car il 
existe des complexes à 12, 13, .....18 électrons et même des complexes riches à 19 et 20 
électrons (dans certains cas rares). 

La richesse ou la pauvreté électronique d’un complexe explique sa réactivité. 

Pour les métaux de transition les plus à droite, notamment ceux de la colonne du cuivre, et 


pour des raisons configurationnelles, la tendance est souvent d'avoir un NEV égal à 14 dans les 
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complexes organométalliques. Ainsi, un grand nombre de complexes à 14 électrons existe +et 
dont la structure est de type MLX (M = Cu, Ag, Au ; L = phosphine par exemple, X = halogène 


par exemple). 


1.5.2. Les équilibres en solution des complexes de métaux de transition 
De manière générale, l’activité d’un complexe métallique, MaL:1Hh, en solution et à une 


température donnée (T), est reliée aux activités des espèces M, L et H par l’équation (1.1) : 
M,LH,)= Bon M)"(L) (H) (LL) 


Par soucis de clarté, les charges des différents ions intervenant dans l’équation ont été omises. 


Best la constante thermodynamique de stabilité de la réaction (1.2) : 


ml 


mM +IL+hH M,LH, (1.2) 


Où m, l et h correspondent respectivement aux coefficients stœchiométriques du métal, du ligand 


et du proton. 
Le logarithme de la constante d’équilibre 8°, est directement relié à l’enthalpie libre AG°de la 


réaction par la relation (1.3) : 


AG° =-RTInf£?, (13) 


La constante globale de stabilité associée à l’équilibre (1.4) peut aussi être exprimée en fonction 


de la concentration des espèces en solution : 
VM,LH, (L4) 


m 


Où les coefficients d'activité Yy, 1,4, ; Ym-71 et 74 rendent compte du comportement non idéal 

des solutions. Le rapport des coefficients d’activité peut être maintenu constant au cours du 

titrage par l’ajout d’un électrolyte support en quantité au moins dix fois supérieure à celle de 

chaque espèce présente en solution. Dans ce cas, la constante £°,, se réduit à une constante de 

stabilité, dite conditionnelle, définie par l’équation (1.5). 
[M,,LH,] 

[m]"ILT (#7 


Lin = (LS) 
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Les constantes de protonation peuvent être exprimées sous la forme de constantes successives 
Koin, définies par l'équation (1.7) et reliées aux constantes de protonation globales par la relation 


(L8). 


LH +H* = LH" (L6) 
LH"? 
Koin = = (L7) 
LEA ]CA*] 
Loin = I] Ki (L8) 
i=1 


L5.3.Description de la structure électronique d’un complexe de métal de transition : 


Pour décrire la structure électronique d’un complexe de métal de transition, il faut combiner 
les orbitales atomiques du métal et les orbitales moléculaires des ligands. A priori, ceci est très 
difficile car il y’a beaucoup d’orbitales. Aussi, les ligands ne sont pas tous équivalents 
électroniquement; ils peuvent contenir des O.M de type sigma et des O.M de type x. Le 


diagramme suivant illustre les interactions de type sigma : 


n OM antiliantes 
o*M-L 


AE (HOMO-LUMO) 


30: , 
A4: (9-n) OM noù liantes 
d À 
18 électrons \ n OM de 
\ 4 type a 
dl 
n OM liantes o M-L 
M L, 

(d'* 3) ML, 


n<9 


Figure 1.13 : Diagramme d'interaction général pour un complexe ML, à 18 électrons 
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On ne peut pas limiter la description électronique des ligands aux seuls effets ©. Les ligands 
possèdent des orbitales moléculaires de symétrie x qui peuvent être vides ou pleines (c’est 
souvent le cas en chimie organométallique). On distinguera alors deux types de ligands [291] : 

Ÿ”__ Les ligands x-donneurs (2 électrons) qui déstabilisent les orbitales d non liantes. 

Y Les ligands z-accepteurs (orbitale moléculaire vide) qui stabilisent les orbitales d non 
liantes. Il est en effet admis que de tels ligands, de par la diversité de leurs 
hétéroatomes(N, O, ….), coordinent différents types de cations métalliques, allant des 
métaux alcalins aux éléments de transition et autres métaux lourds, et s'utilisent dans des 


domaines aussi divers et disparates que la catalyse, la biologie, la corrosion, … [30]. 


L.6. Les complexes de coordination en biologie 


Selon Brown et al. 1980 [31], de nombreux médicaments possèdent des propriétés 
pharmacologiques et toxicologiques modifiées lorsqu'ils sont administrés sous forme de 
complexes métalliques. Le cation le plus largement étudié à cet égard est probablement le Cu°* 
car de nombreux complexes de cuivre de faible poids moléculaire se sont révélés bénéfiques 
contre plusieurs maladies telles que la tuberculose, les rhumatoïdes, les ulcères gastriques et les 
cancers. Vingt acides aminés naturels constituent les éléments constitutifs des protéines, qui sont 
des espèces chimiques indispensables à de nombreuses fonctions biologiques, comme en 
témoigne le rôle des enzymes. 

Les métaux impliqués dans des processus biologiques sont, soit liés directement aux 
protéines, soit à des cofacteurs (comme les porphyrines ou les cobalamines), soit à des clusters, 
eux-mêmes ligands de protéines. Les principaux atomes ligands des métaux sont l’azote, 
l’oxygène, le soufre et le carbone [321]. 

Les complexes de cuivre à base d’azote-donneur sont de grands intérêts car ils peuvent imiter 
les sites actifs de nombreuses enzymes contenant du cuivre et aussi pour leur utilisation 
potentielle en tant que métallo-médicaments. Le cuivre en tant qu'élément essentiel pour 
l'homme, de par son activité bio-essentielle et sa nature oxydante, a amené de nombreux 
chimistes inorganiques à aborder les complexes de Cu (Il) dans des applications médicales. Il 
existe plusieurs articles scientifiques dans la littérature concernant les activités biologiques de 
ces complexes synthétiques. Les complexes de cuivre (Il) peuvent interagir avec les acides 
nucléiques par intercalation ; toutefois, leurs interactions dépendent de leur structure 
moléculaire. À cet égard, la conception de biomatériaux fonctionnels a fait l'objet d'une attention 


particulière en raison de leurs applications potentielles en tant que matériaux pharmaceutiques 
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[33]. Bien que le zinc, qui possède un orbital d complet, ne soit pas au sens strict un métal de 


transition, 1l partage beaucoup de propriétés bio-inorganiques avec les métaux de transition. 
> Les complexes du cobalt 


Les complexes du cobalt ont montré des effets anti-paludéens où ils visent le plasmodium 
parasite (Figure 1.14-a). Ces études étaient in vitro et ne montrent pas la toxicité au parasite 
l'hôte. D'autres études ont montré des complexes du cobalt pour avoir activité anti-bactérienne 
(Figure 1.14-b). Ces résultats ont aussi été faits in vitro. Le cobalt est un redox métal actif oùles 


+2 et +3 états de l'oxydation sont très communs. 


à N2H / 
Z ÿ H 
NH 
2 NO 
R=NTO TN Nr 
N Na 
TT —Q CI N 
HN 
DS à " 7 
(o) H2N 
(a) | d 


Figure 1.14 : Exemple de complexes de cobalt utilisés dans les études antipaludiques (a) et 


antibactériennes (b) 


> Les complexes du zinc 
Les complexes du zinc (Figure 1.15) ont été montrés pour lier fortement à un récepteur HIV 
CXCRA4, utilisé pourla fusion de la membrane et d'où inhibe l'étendue du virus. Les complexes 
ont été montréspour être dix fois plus actifs que le ligand libre. D'autres complexes du zinc ont 
été montrés pour avoir l'activité anti-paludéenne et aussi activité anti-bactérienne. 
o [e) 
S à NS NN 
| L AT u ji di 
22 \ Rs DA H 22 \ ne DZ H 
PON| Zn Zn |N7 | POSN| Zn Zn |N7 | 
K N D NA K ! N' EN 


R NE 


1R=n-Bu 6R= n-Bu 


2R=1-Bu TR= TC 
3R=n-octyl F 
4R=2-hydroxyphenyl 

SR=2-azidoethyl 
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Figure 1.15 : Complexe du ZnDPAantiplasmodiquesynthétisé par Rice[34] 
> Complexe de palladium 


Le complexe de palladium (Figure 1.16) a montré l'activité anticancéreuse semblable à leur 
équivalents du platine et dans quelques cas a été plus actif que cis platin. Pd (ID) est un d8 ion qui 


a la géométrie du plan carrée autour du métal, semblable à Pt (ID. 


NH, 


ee 


/ NT 


CI 7 CI 


Figure 1.16 : Exemple d'un complexe de palladium utilisé dans une étude anti-cancer 


L.6.1. Affinité des métaux pour les protéines 


Tout au long de la série des métaux de transition (de Mn°* à Cu?*), la taille des ions diminue, 
tandis que l’affinité de l’ion métallique pour les électrons augmente (ce qui est inversé pour le 
zinc). En conséquence, l’affinité relative du métal pour les atomes non-métalliques ou les anions, 
qui peuvent donner ou partager des électrons, augmente. Ainsi, la série d’Irving-Williams [35] 
classe les métaux en fonction de leur capacité d’interaction avec les ligands [36]. Ca?*< Mg?*< 
Mn’*< Fe*< Co?*< Ni2* < Cu?*> Zn°?* 

En général, l’affinité des métaux pour les métalloprotéines suit l’ordre défini par cette série. 
Ainsi, le cuivre et le zinc ont souvent des affinités plus fortes pour la majorité des 
métalloprotéines que les autres métaux. Ces ions doivent donc être tenus éloignés des sites 
réservés aux cations de plus faible affinité. Des systèmes de transport et de stockage permettent 


de maintenir une concentration extrêmement faible de ces ions sous forme libre dans le 


cytoplasme des cellules [36]. 


1.6.2. Protéines à cuivre 


ES 


Les protéines à cuivre sont des enzymes qui possèdent des sites actifs présentant un ou 
plusieurs atomes de cuivre [37]. Ce dernier a un rôle de cofacteur et il intervient donc dans de 


nombreux processus tels que le transfert d’électrons, le transport de l’oxygène, le métabolisme 
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du fer, la régulation des neurotransmetteurs et des neuropeptides, la synthèse des phospholipides, 
le fonctionnement du cerveau, la pigmentation, etc. 

L'utilisation du cuivre par le vivant est étroitement liée aux propriétés d’oxydo-réduction du 
couple Cu*/Cu?“*dont le potentiel standard redox, est de 160 mV. Le cuivre a aussi une capacité à 
former des liaisons de coordination stables avec les groupements thiolate et imidazole. Les 
protéines à cuivre peuvent se classer en plusieurs types selon la structure dusite liant le cuivre. 
Certaines enzymes sont constituées de sites plus complexes de type, CuaA( Figure 1.17) oùles 
atomes de cuivre sont liés par des histidines, une méthionine et un carbonyle d’un acide aminé et 
sont pontés par un atome de soufre d’une cystéine, ou de type Cugoù le cuivre est coordonné par 
trois histidines dans une géométrie trigonale (Figure 1.17) ou plus complexe encore, avec 
d’autres ions métalliques tels que le fer, le manganèse ou le zinc. Enfin, le site Cuz que l’on 
rencontre dans l’enzyme oxyde nitreux réductase est quant à lui constitué de quatre atomes de 


cuivre, entourés par 7 histidines et reliés par un atome de soufre pontant (Figure 1.17)[38]. 


Glu Hi 
Cys  N His À His 
. : TT 
! His | Ms HO—Cu 
His ” F Ho,,, His 
nd cu—© “ : CUBE: Cu”, . 
LR , His\ His 
Met s PAS 4 | 
Cys ‘Fe HisŸ Vue 


Figure 1.17 : Centre de type Cua, CugetCuz 


L.6.3. Propriétés antibactériennes des métaux de transition 


Les ions métalliques ont été explorés en médecine pendant longtemps. Selon la littérature, la 
vaste application des métaux en médecine comprend des composés antimicrobiens et 
anticancéreux, des agents antithritiques, des inhibiteurs enzymatiques, etc. [391]. 

Depuis l’antiquité, le cuivre est utilisé par l’homme notamment pour soigner les infections et 
prévenir les maladies. Avant la découverte des micro-organismes, l’homme a utilisé des 
préparations à base de cuivre pour les maux de gorge, éruptions cutanées et pour l’hygiène 
quotidienne. Au XIX° siècle, après la découverte du lien de causalité entre le développement de 
germes pathogènes et la déclaration des maladies, de nombreux scientifiques se sont intéressés à 
l’exploitation des propriétés antibactériennes du cuivre. Actuellement, le cuivre est utilisé par 
l’industrie pharmaceutique, dans des applications allant des antiseptiques et antimycosiques aux 


produits de soins et d’hygiène (crèmes, ampoules d’oligo-éléments, etc.). 
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En 2015, des complexes de Co (I), Ni (ID, Cu (I) et Pd (I) de 4-amino-S pyrimidine 
carbonitrile (APC) ont été synthétisés et caractérisés par Mohamed et al. [39]. La Figure 1.18- 
aprésente le mode d’interaction d’APC avec Co (ID) et l’efficacité antimicrobienne d’APC et de 
leurs complexes. L’efficacité antimicrobienne de ces complexes pourrait être liée à l’inhibition 
de plusieurs enzymes structurelles joue un rôle clé dans la vie voies métaboliques des 
microorganismes, le ligand APC et ses complexes métalliques ont montré une activité modérée 


contre les souches bactériennes par rapport au médicament standard (tétracycline)(Figure 1.18- 
b). 


Inhibition Zone (mm) 


sc 


APC Co(i) Ni(Il) Cut) Pd(it) Tetracycline 


a b 


Figure 1.18 : (a) Mode d’interaction d’APC avec Co (ID), (b) Efficacité antimicrobienne 


d’APC et de leurs complexes 


1.7. Conclusion 


Dans ce chapitre, nous avons rapporté des généralités sur les complexes de coordination. 
Nous avons, également, fait le point sur les différents types des ligands et leurs méthodes de 
coordination avec les ions métalliques ainsi que leurs réactivités. Les ions métalliques de 
coordination ont reçu une attention considérable pour les chimistes en raison de leurs utilisations 


en synthèse de complexes métallo-organiques d’intérêt biologique. 
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Chapitre II 


Sulfonamides, Oxazolidinones et N-sulfonylimines 


IL.1. Introduction 


Dans ce chapitre, nous présentons des généralités sur les sulfonamides et leurs dérivés N- 
sulfamoyl oxazolidinones et N-sulfonylimines, en donnant leur mode de coordination au sein 
des complexes métalliques et organométalliques, et l’importance de cette interaction sur les 


activités biologiques. 
IL2. Généralités sur les sulfonamides et leurs complexes de coordination 


Un sulfonamide est un membre d'une classe de composés chimiques contenant le groupe - 
SO2N-. La classe comprend plusieurs groupes de médicaments utilisés dans le traitement des 
infections bactériennes, du diabète sucré, de l’œdème et de l’hypertension [1]. 

Les sulfamides sont des analogues à l’acide p- aminobenzoïque (PABA). Le prontosil, un 
colorant azoïque métabolisé in vivo en sulfanilamide, a été le premier sulfonamide d'importance 
clinique. De nombreux sulfamides ont été synthétisés, mais leur activité antimicrobienne ne 
diffère que légèrement, tandis que leurs propriétés pharmacocinétiques varient. Les sulfamides 
ont été classés en fonction de leurs taux d’excrétion en activités : courtes, moyennes, 
intermédiaires ou longues et ultra-longues. Les sulfamides à courte durée d’action sont excrétés 
dans l’urine de fortes concentrations et ont donc été particulièrement utiles dans le traitement des 
infections des voies urinaires. Parmi les sulfamides à action brève les plus couramment utilisés, 
la sulfadiazine a une faible solubilité dans l'urine alors que la sulfamérazine et la sulfaméthazine 
et leurs conjugués acétyle sont très solubles. 

Les sulfamides sont des antibiotiques bactériostatiques synthétiques à large spectre contre la 
plupart des organismes gram-positifs et de nombreux organismes gram-négatifs. Cependant, de 


nombreuses souches d'une espèce donnée peuvent être résistantes. Les sulfamides inhibent la 
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multiplication des bactéries en agissant comme inhibiteurs compétitifs de l'acide p- 
aminobenzoïque dans le cycle du métabolisme de l'acide folique. La sensibilité bactérienne est la 
même pour les différents sulfamides et la résistance à un sulfamide indique une résistance à tous. 

Ces médicaments sont largement répondus en médecine et ont été les premiers agents utilisés 
dans le traitement de l'infection bactérienne pour les maladies humaines et animales. 
Aujourd'hui, l'utilisation de sulfamides est limitée au traitement de maladies spécifiques en 
médecine humaine, telles que les infections des voies urinaires. Cependant, les sulfamides sont 
plus souvent rencontrés en médecine animale. La présence de certains résidus dans les produits 
animaux présente un risque potentiel pour la santé en raison de leurs propriétés allergiques [21]. 
De plus, certaines personnes présentent une hypersensibilité aux résidus de médicaments et / ou 
aux faibles concentrations des résidus de médicaments pouvant produire des bactéries 
génétiquement modifiées résistantes au traitement médicamenteux existant [3]. En outre, une 
étude du Centre national de recherche toxicologique a indiqué que la sulfaméthazine pourrait être 


un cancérogène de la thyroïde [4]. 
IL.2.1. Historique 


Les sulfamides ont été synthétisés pour la première fois par Gelmo et al.[5] en 1908 alors 
qu'ils effectuaient des recherches sur les colorants azoïques. Juste après ces travaux, Hoerlien et 
al. [6] ont découvert des colorants contenant le groupe sulfanyle ayant une affinité pour les 
protéines de la soie et de la laine. 

Cependant, jusqu'en 1932, les propriétés thérapeutiques des sulfonamides n'étaient pas 
reconnues. Les scientifiques allemands Domagk et al.[7] ont découvert que le prontosil avait une 
activité antibactérienne in vivo prononcée. Ils ont observé que les souries atteintes de septicémie 
à streptocoque pouvaient être guéries avec du prontosil. Domagk a également découvert que le 
prontosil était rapidement réduit dans la cellule en sulfanilamide et qu’il s’agissait en fait du 
sulfanilamide et non du prontosil (Figure 2.1), qui était le véritable antibiotique [8].Cette 
nouvelle classe de sulfamides était la seule antibactérienne efficace jusqu’à ce que la pénicilline 
soit disponible au début des années 1940. Malgré les avantages des sulfamides, ils se sont révélés 


inefficaces contre des infections telles que la salmonelle (organisme responsable de la typhoïde) 


[9]. 
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NH; 
Î 
HN N=N = —NH, Réduction 
Prontosil l 
NH; 
ï 
HN NH; + HN ac 
O 


Sufanilamide 
Figure 2.1 : Sulfanilamide de prontosil 


De nombreux sulfamides différents ont été synthétisés à la fin des années trente. On a 
découvert qu'un grand nombre d'entre eux possédaient une activité antibactérienne considérable 
pour une variété de bactéries streptocoques et pneu monocoques. Plusieurs sulfapyrimidines 
introduites en 1941 [9] se sont avérés posséder une activité antibactérienne considérable et une 
toxicité inférieure à celles des sulfamides antérieurs. Cette avancée a conduit à la synthèse de 
nombreux nouveaux sulfamides. Aujourd’hui, il existe plus de 5000 sulfamides, mais seulement 
33 d’entre eux ont été introduit pour un usage médical général. 

L'utilisation des sulfamides en tant que médicaments remonte au début du XX® siècle, lorsque 
la découverte de l’usage médicinal des sulfamides et de leurs dérivés a marqué un tournant dans 
l’histoire de la chimiothérapie (Figure 2.2). Il s'agissait de la première enquête de synthèse des 
molécules organiques en tant que médicaments potentiels pour combattre les infections 
véhiculées par le sang. 

L'importance des sulfonamides dans l’industrie pharmaceutique a même été mise en évidence 
lors de la recherche dans la base de données lorsque le composé 5-acétamido-Il, 3, 4-thiadiazole- 
2-sulfonamide était localisé. L’acétazolamide, ou Diamox comme il est communément appelé, 
est utilisé dans le traitement de nombreux types de glaucome (pression élevée dans la chambre 
antérieure de l’œil ; la chambre antérieure est la partie la plus avancée de l'œil, derrière la 
cornée). Il est également utilisé comme diurétique dans certains problèmes tels que l'insuffisance 
cardiaque ou lorsqu'un liquide est retenu dans le corps. Ce médicament empêche l'action d'une 
enzyme, l'anhydrase carbonique, dans tout le corps. Dans l'œil, cela entraîne une diminution de la 
production de fluide sécrétée dans la chambre antérieure de l'œil, ce qui entraîne une diminution 
de la pression. Il agit dans les reins pour augmenter la quantité de bicarbonate excrété dans 
l'urine et augmente ainsi l'élimination du sel et de l'eau du corps. La façon dont cela fonctionne 


sur le système nerveux central n'est pas encore bien comprise [9,10]. La relative insolubilité de la 
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plupart des sulfonamides, en particulier de leurs métabolites acétylés, peut provoquer une 
précipitation dans les tubules rénaux. Les analogues les plus solubles, tels que le sulfisoxazole et 
le sulfaméthoxazole, doivent généralement être choisis et le patient doit être bien hydraté. Pour 
éviter la cristallerie et les lésions rénales, la consommation de liquide devrait être suffisante pour 


produire un débit urinaire de 1200 à 1500 m1 / jour [9,101]. 


Figure 2.2 : Interactions 2D d’un sulfonamide avec les acides aminés du site actif de 


pdb :1AJ0 (CIRCULARLY PERMUTED (1-3,1-4)-BETA-D-GLUCAN 4-GLUCANOHYDROLASE CPA16M- 
127) 


IL2.2. Travaux antérieurs sur les médicaments à base des complexes de coordination 


contenant des sulfonamide 


La présence de plusieurs sites donneurs potentiels, par exemple l'azote sulfonamide, un azote 
amino, deux atomes d'azote pyrimido et deux atomes d'oxygène sulfonyle en font un agent 
complexant polyvalent avec des ions métalliques. Les dérivés de sulfonamide et leurs complexes 
métalliques possèdent une activité antibactérienne. Quelques complexes métalliques des 
sulfonamides sont encore utilisés dans la médecine topique comme antiseptiques, bien qu'il soit 
connu que les interactions d'ions de métaux lourds avec des biomolécules peuvent avoir des 
effets potentiellement toxiques. Les complexes métalliques des médicaments à base de 
sulfonamide ont acquis un rôle important dans la chimie de coordination. Cependant, les études 
rapportées sur les complexes métalliques ne sont ni systématiques n1 exhaustives (Figure 2.3) 


[10]. 
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HN 


Figure 2.3 : Intéraction des métaux de coordination M?*avec l'azote sulfonamide 
8 


Selon Behrens et al, 1986[11], on sait que plusieurs chélates métalliques possèdent une 
activité antibactérienne, anti-fongicide, antivirale et anticancéreuse. Dans plusieurs cas, les 
chélates de métal ont été plus antimicrobiens que les agents calcifiants eux-mêmes. On sait 
également que certaines drogues agissent par chélation ou par métalloenzymes inhibitrices, mais 
pour la plupart des médicaments qui agissent comme des ligands potentiels, de nombreuses 
études sont réalisée pour déterminer comment les liaisons des métaux influés sur l’activité des 
médicaments d’autres sont en cours. 

Cependant, les ions métalliques tels que Fe, Co, Ni, Cu, Zn et Cd jouent un rôle important dans 
la chimie bio-inorganique. Ces composés sont extrêmement importants, bien que quelques 
articles aient été publiés sur les complexes de métaux de transition de ligands apparentés [12]. 

Il est à noter que, outre l’intérêt stéréochimique des complexes métalliques de sulfonamide, 
leur utilisation clinique est importante. Par exemple, les complexes d'argent et de zinc ont déjà 
trouvé une application en tant que médicaments topiques pour la guérison des brûlures. 

Les sulfamides fournissent des sites de liaison potentiels pour les ions métalliques et toute 
information sur leurs propriétés de coordination est importante pour comprendre le rôle des ions 
métalliques dans les systèmes biologiques. L'activité biologique prédominante de divers 
complexes sulfonamides métalliques a été rapportée avec des activités antibactériennes [13], 
anti-tumorales, diurétiques,  anti-carboniques  anhydrase [14],  hypoglycémiques, 
antithyroïdiennes, protéases et plusieurs autres activités [15]. Les complexes de sulfanamides Ag 
(D et Zn (I) sont utilisés avec succès depuis plusieurs décennies pour la prophylaxie et le 
traitement des infections de brûlures microbiennes et fongiques [16].Dans la Figure 2.4 
Mastrolorenzo et al.[17]proposons que des ligands agissent de manière monodentée dans les 
nouveaux complexes Ag (1) (A), le système donneur étant constitué par l'azote sulfonamidique 
du fragment p-phénylène-SO2NH-thiadiazole. D'autre part, les mêmes ligands agissent 


probablement de manière bidentée dans le Zn (I) (B). 
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Figure 2.4 : Nouveaux complexes de sulfonamidique avec Ag (1) et Zn (I) 


De plus, il a été constaté que plusieurs complexes de sulfonamide de Cu (I), Ce (D), Bi (D), 
Cd (D) et Hg (I) ont révélé une activité antibactérienne. Les complexes de sulfonamides, ainsi 
que leurs dérivés, trouvent également une application en chimiothérapie. 

Différents dérivés de sulfonamide substitués contenant un cycle oxazolidinone sont une 
classe essentielle de composés bioactifs ayant un large spectre d'activités pharmacologiques, 
parmi lesquels les N-sulfamoyloxazolidinones, qui sont des composés attrayants combinant un 
pharmacophore oxazoldinone et un radical sulfamoyle. Ils se sont avérés efficaces comme agents 
antibactériens et ont été utilisés comme précurseurs dans la synthèse de 2- 
chloroéthylnitrososulfamides CENS. 

Certains auteurs ont modifié ou remplacé le cycle de l'oxazolidinone, ce qui a entraîné une 
amélioration considérable de la puissance des oxazolidinones [18-22]. Tant de critiques 
intéressantes ont été publiés par différents groupes de recherche sur les bactéries antibactériennes 
à base d'oxazolidinone [18,23-26]. 

Ces hétérocycles présentent un intérêt considérable en chimie pharmaceutique et médicinale 
et constituent des classes de composés hétérocycliques avec une large gamme d'activités 


biologiques [27]. En outre, l’importance de l’imidazolidine-2-one a été soulignée dans la 


33 


Chapitre II Sulfonamides, Oxazolidinones et N-sulfonylimines 


synthèse organique. Les oxazolidines constituent une classe très importante de composés 


hétérocycliques présentant des activités biologiques intéressantes [28]. 


IL.3. Les oxazolidinones 


L'étude des systèmes hétérocycliques à cinq chaînons a connu un développement 
considérable, dû à la mise en évidence de leurs activités biologiques. Les oxazolidinones sont 
parmi une large variété des hétérocycles étudiés pour développer de nouvelles molécules 
bioactives. L’oxazolidinone (Figure 2.5) est une structure la plus largement utilisée en synthèse 
asymétrique et qui se retrouve un peu moins fréquemment dans les produits naturels. Ce 
composé a été introduit pour la première fois par Evans en 1981[29], cette molécule présente une 


activité remarquable in vitro et in vivo. 


R; 


R KR 


OL > SHN 


O 


Figure 2.5 : Structure générale d’oxazolidinone 


L'importance des oxazolidinones se situe au niveau du rôle qu’elles jouent comme auxiliaires 
chiraux. En effet, elles représentent un groupe important d’auxiliaires utilisés pour contrôler la 
stéréosélectivité de nombreuses transformations chimiques diversifiées comme; les réactions 
d’aldolisation, les réactions de cycloaddition de Diels-Alder, et les réactions d’alkylation des 
énolates [30]. Parmi les auxiliaires de type oxazolidinone les plus connus, on trouve les 
auxiliaires d’Evans (Figure 2.6) [31]. Le groupe de Sibi a développé un nouvel auxiliaire 


particulièrement efficace pour les réactions radicalaires [32]. 
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Figure 2.6 : Exemple des auxiliaires d’oxazolidinones les plus connus. 


Le N-aryloxazolidinone est un motif structural important pour l'utilisation médicinale. 
Plusieurs médicaments cliniques portant ce fragment ont été développés comme l’éperzolide 
et le linézolide.….etc. J. Li et al. [33], ont décrit une réaction d’arylation des oxazolidinones avec 
le bromure d'aryle en présence de K2CO3 comme base pour former les Naryloxazolidinones. Ils 
ont utilisé le N, N-diméthylglycine comme un excellent ligandpour le couplage catalysé par du 


Cul à 120 °C après 24 h (Figure 2.7). 


O 
— Vo 5 mol% Cul, 
B EL % DMF 
Tr + 
——. l HN 7 K2CO;, DMF_ 
R;120 °C 


R; 


R;=H, Ph, i-Pr R;=H, CH;-OH  X=H, MeOC, CN, NO;, SMe, Me, OMe 
Figure 2.7 : Synthèse d’aryl-oxazolidinone 


IL.3.1. Activités pharmacologiques des oxazolidinones 


Les oxazolidinones sont largement utilisés dans la synthèse des produits pharmaceutiques, les 
pesticides et les cosmétiques [34]. La classe des oxazolidinones est la première classe 
d'antibiotiques qui soit apparue au cours de ces 30 dernières années. Cette classe de composés 
montre une excellente activité contre les bactéries à Gram positives, aussi bien que plusieurs 
anaérobies et Mycobactérie tuberculoses[35]. Dans un premier temps, elles sont utilisées comme 
antidépresseurs inhibiteurs de la monoamine oxydase (IMAO) pour le traitement de dépression 


et de la maladie de Parkinson, puis comme myorelaxant, anticonvulsivant et tranquillisant et 
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enfin pour leur activité fongicide vis-à-vis de la tomate, nous citons à titre d’exemple le 


toloxatone un inhibiteur de la monoamine oxydase (IMAO) (Figure 2.8). 


HO 


O 
Figure 2.8 : Structure chimique de toloxantone. 


À la fin des années 80 ; les oxazolidinones ont été utilisées pour certains agents pathogènes 
humains. En 1987 Slee et al. [36], décrivent pour la première fois des composés appartenant à la 
classe des oxazolidinones, nommés DuP 105 et DuP 721 (DuPont de Neumours), développés 


afin d'être actifs envers les germes à Gram positif (Figure 2.9) 


O 
O Vo Yo H 
Mel Ÿ VX À CAS 
O O 
Pa Me PE 
Dup 105 Dup 721 


Figure 2.9 : Structures chimiques de Dup 105 et Dup 721. 


Le linézolide [37,38] est le premier et le seul médicament à base d’oxazolidinone qui a 
étéapprouvé par la FDA en 2000 pour le traitement de la polychimio-thérapieinfections 
bactériennes à Gram positif résistantes. Dans un court laps de temps deson lancement, des 
souches de S. aureus et d’Enterococcus résistantes au linézolide faecium{[39]ont commencé à 
émerger. Il y a un besoin urgent d’antibiotiques plus récents et améliorés à large spectre 
antimicrobien qui sont également actifs contre les souches résistantes au linézolide. De 
nombreux efforts ont été déployés pour améliorer la puissance et le spectre antimicrobien de ce 


médicament [40-42] (Figure 2.10). 
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Linézolide Eperozolide 
Figure 2.10 :Structures du linézolide et de l’épèrozolide 


IL.3.3. Complexes de coordination avec l'oxazolidinones 


En plus d'être des composés antimicrobiens, les dérivés d’oxazolidinones ont des ligands 
potentiels à partir desquels peuvent êtres générer des complexes métalliques. Les complexes 
métalliques d'oxazolidinones étaient connus et utilisés dans la catalyse à savoir. la réaction de 
Diels-Alder, [43-45]Jréaction aldol, [46-48] réaction de Henry, [491 réaction de Mukaiyama — 
Michael, [50] diamination des alcènes [51-S2]ou pour générer des MOFT[53]principalement en 
fonction du choix du métal et de la coordination environnement rendu. Bien que les complexes 
du métal à base d’oxazolidinone ont été utilisés en catalyse et hautement d'intérêt thérapeutique à 


savoir agents antimicrobiens ou anticancéreux. 


O 


JS 


Figure 2.11 : Synthèse des complexes métalliques à base d’oxazolidinone 1-3[54] 
IL.4. Les bases de Schiff et les N-sulfonylimines 


Les sulfonamides et leurs composés dérivés de la base de Schiff sont largement utilisés pour 
antibactérien[55], anti-tumeur[56], diurétique[57], anti-anhydrase carbonique (anti-CA)[58], 
hypoglycémiant[59], anti-thyroïdien[60] et inhibiteur de protéase[61] activités. 
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De nombreux médicaments possèdent des propriétés pharmacologiques et toxicologiques 
modifiées, lorsqu'ils sont administrés sous la forme de leurs complexes métalliques. Les ions 
métalliques les plus étudiés à cet égard sont le cobalt (ID), le cuivre (I), le nickel (ID) et le zinc 


(ID [62,63]. 
II.4.1. Les bases de Schiff 


Une base de Schiff est toute molécule comportant une fonction imine (C=N). Le chimiste 
allemand Hugo Schiff fut le premier à avoir synthétisé ce type de composés en 1864; c’est 
ainsi qu’ils gardent son nom jusqu’à présent (Bases de Schiff) [64]. Ils sont obtenus 


par condensation d'une amine primaire avec un aldéhyde ou une cétone (Figure 2.12). 


OH 
H R; 
Ne + R—NH, H—C—R + H—C—R + HO 
l NHR NR 


Figure 2.12 : Processus réactionnel d’obtention d’une base de Schifff65] 


IL.4.2. Les N-sulfonylimines 


Plusieurs réactions de condensation des sulfonamides avec les aldéhydes pour accéder aux N- 
sulfonylimines ont été décrites dans la littérature [66-68]. Durant les dernières décennies, les N- 
sulfonylimines ont gagné une importance, particulière pour les chimistes organiciens, et ce, dans 
les domaines de la méthodologie, de l’industrie ainsi qu'en synthèse totale de molécules 
naturelles. Plus encore, leur utilisation en synthèse asymétrique constitue en soi un défi pour les 
groupes de synthèse organique. La N-sulfonylimine (Figure 2.13) est une classe spéciale d’imine 
contenant le motif sulfonyle [69], utilisé comme précurseur en synthèse de composés 
biologiquement actifs [70]. En effet, il est à noter que le groupement sulfonyle accroît 


l’électrophilie de la liaison C=N ce qui permet d’atteindre d’excellents rendements avec 


l’addition d’organométallique. 
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R102 N 


Figure 2.13 : Structure générale de N-sulfonylimine 


La double liaison carbone-azote C=N du N-sulfonylimine a beaucoup attiré l'attention des 
chimistes comme un excellent accepteur de radicaux [71], et d'excellent précurseur pour 
lapréparation des aziridines[721, des imidazolines [73], des N-sulfonyloxaziridines [74]qui ont 


une utilité d’oxydants chiraux et dans la synthèse asymétrique des dérivés de f-aminoacide[751. 


IL.4.3. Propriétés complexantes des bases de schiffet des N-sulfonylimines 
IL.4.3.1. Coordination à l’état solide 


La découverte des propriétés complexantes de certaines bases de Schiff par V. P. Pfeifferen 
1931 a conféré à ces dernières une importance particulière et sans cesse grandissante, tel que le 
prouve la multitude de travaux dont elles font l’objet au fil du temps [76]. 

Les données de la littérature ont permis à titre d’exemple d’établir des complexes 
polymères de type ML: avec les métaux de la première série de transition. Selon la nature du 
milieu et sa température, ces complexes s’avèrent adopter différentes géométries, octaédriques, 
plan carrées ou même un équilibre entre les deux [77-81]. 

En effet, V. P. Pfeiffer et T. Tsumaki ont été les premiers à élaborer des complexes de métaux 
de transition avec des ligands tétradentés à pont éthylidène, usuellement appelés Salen (Figure 
2.14) [82-85]. Il s’agit là du début d’une méthode de synthèse des complexes tétradentés, 


généralisée ultérieurement par les mêmes auteurs. 


À 
N N— —= = 
M(OA“h 
DEL 
OH HO 
M=Co(Il), Ni(ID), Cu(ID) et Zn(Il) 


Figure 2.14 : Synthèse des premiers complexes dérivés de base de Schiff de type Salen 


Les complexes de métaux de transition avec des sulfonamides de dérivé bases de Schiff 


donneurs d'oxygène et d'azote présentent un intérêt particulier en raison de leur capacité à 
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posséder des propriétés géométrique inhabituelles [86]. La Figure 2.15 présente une structure 


proposée des complexes métalliques à base d’un sulfonamide de dérivé bases de Schiff. 


D.” TT — 
PHÉRRS ST 
—N N 
Le H N AT H — 


M=Cu(Il); x=0; n —2 
Zn(I); x=0; n =3 
Cdi); x=0; n =5 


Figure 2.15 : Structure proposée de [M(N,N0-bis(5-(4-sulfanilamidophenylazosalicylidene)- 
ethylenediamine)] synthétisé parA.-N.M.A. Alaghaz et al [87]. 


En effet, les N-sulfonylimines constituent les premiers systèmes ayant fait l’objet d’études de 
coordination vis-à-vis de cations métalliques divers tels, le nickel(ID), cuivre(I), fer(H). 
La structure de ces ligands s’avère d’autant plus intéressante qu’elle renferme, outre 
N-sulfonylimines, un groupement C=N, augmentant leur basicité et leurpermettant de 
coordonner, avec des constantes de stabilité élevées, différentes entitésmétalliques, allant des 


métaux de la première série de transition aux métaux lourds, tels que l’uranyle ou le tungstène. 


IL.4.3.2. Les complexes de sulfonamides et leurs composés dérivés de la base de Schiff en 


biologie 


Certains complexes métalliques sulfonamides ont beaucoup attiré l’attention du fait que leur 
activité biologique était supérieure à celle des ligands sulfonamides ou les sels métalliques [88]. 

Les ligands imine et imine réduite et leurs complexes métalliques sont également utilisés dans 
divers domaines d’application en chimie des polymères, en raison de leurs propriétés 
magnétiques, de leurs propriétés optiques, comportement thermique, chimie médicale, chimie 
des catalyseurs...etc. [89-100]. 
En outre, les ligands de sulfonamide de dérivé de base de Schiff et leur métal correspondant 
complexes se sont énormément étendus et comprennent un vaste domaine de composés 
organométalliques et divers aspects de la chimie bio-inorganique [55-60]. Des ligands de la base 
de Schiff ont été rapportés pour montrer une variété d'actions biologiques. 

Depuis quelques années, on assiste à travers la littérature scientifique à un intérêt 
concernant les bases de Schiff qui présentent des activités biologiques spécifiques en 
général reliées à la présence de ce groupement fonctionnel. 
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La lutte contre les maladies infectieuses constitue, de nos jours, un enjeu planétaire de 
la plus haute importance. Les problèmes infectieux demeurent en effet des plus préoccupants, 
du fait du surcout hospitalier qui en découle, de l’émergence de bactéries multi-résistantes 
(BMR) et du taux de mortalité dont elles sont responsables. Pour cela, la recherche de 
nouvelles molécules dotées de propriétés antibactériennes, antifongiques et anti-inflammatoires 
demeure un sujet d’actualité, qui revêt une importance sans cesse grandissante avec le temps. 
Les composés organiques de type sulfonamide de dérivé de bases de Schiff et leurs complexes 
métalliques constituent de nos jours des systèmes dont l'importance et l'efficacité thérapeutique 
en tant qu’agents anti-viraux et anti-cancéreux demeurent avérée. 

Nous citerons dans ce qui suit, quelques-unes des références bibliographiques les plus 
récentes, ayant porté sur l'élaboration de bases de Schiff analogues aux nôtres et sur la mise en 
évidence de leurs propriétés bioactives ainsi que celles de leurs complexes avec différents 
cations métalliques. 


c 


o Gly122(A) 


Meté1(A) — Fi 


Cys163(A) 


- -Tyr204(B) 
Phe256(B}) 
Hisi21(A) 


-_Trp206(B) 


IGFW_xray 
Figure 2.16 : Interactions de l'inhibiteur à base d'isatine MSI présent dans la structure 
cristalline 1GFW. Les liaisons hydrogène sont représentées par des lignes vertes. Les contacts 
hydrophobes sont représentés par des hémisphères radiaux. Les interactions de liaison ont été 


calculées avec le programme LigPlot.20 [1011]. 
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ILS. Conclusion 
Selon les travaux présentés dans cette partie, il en ressort que l’insertion d’un ligand dans la 
structure d’un complexe métallique ou organométallique contribue à l’amélioration de l’activité 


biologique et constitue une stratégie porteuse dans cet axe de travail. 
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Chapitre III 


Etudes théoriques 


IIL.1. Introduction 


La modélisation moléculaire est un outil destiné aux chercheurs préoccupés par la structure et 
la réactivité des molécules. La connaissance de la structure des édifices moléculaires permet de 
comprendre le résultat d’une transformation physique, chimique ou biologique ou encore de 
prévoir de telles transformations. 

Dans ce chapitre nous avons présenté les méthodes théoriques DFT permettant l’étude des 
géométries et de la réactivité de nos molécules, et aussi le principe général du Docking 
moléculaire pour interpréter les interactions intermoléculaires Ligands-récepteur protéique par 


docking moléculaire. 
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IIL.2. Modélisation moléculaire 


La recherche et la synthèse de nouveaux composés chimiques sont aujourd’hui souvent 
associées à une étude par modélisation moléculaire. La modélisation moléculaire est une 
technique permettant, non seulement de représenter les propriétés et les réactions chimiques mais 
aussi de manipuler les modèles des structures en deux ou trois dimensions. 

La modélisation moléculaire implique l'utilisation des méthodes de calcul théoriques 
(mécanique moléculaire, dynamique moléculaire, mécanique quantique ab-initio ou semi- 
empirique...) permettant de déterminer la représentation graphique en 2D ou 3D de la géométrie 
ou de la configuration des atomes d'une molécule, et d'évaluer les propriétés physico-chimiques 
de la molécule étudiée. 

Ces approches présentent certains avantages par rapport à leurs homologues 
expérimentaux, telle que la compréhension des caractéristiques des contacts intermoléculaires, 
contrairement aux mesures expérimentales. Dans le cadre de cette thèse, ces outils théoriques 
sont employés pour deux raisons : 

> La première consiste à mimer les molécules synthétisées, par les méthodes de la théorie de 

la fonctionnelle de la densité (DFT), et pouvoir par la suite étudier leurs 
comportements en déduisant ainsi leurs différentes propriétés physico-chimiques (pour 
les comparer avec les résultats expérimentaux). 

> La deuxième raison est destinée à prévoir un système idéal via un criblage virtuel 

présentant les différents contacts potentiels entre les inhibiteurs et leurs récepteurs 


protéiques (pour expliquer les résultats enzymatiques). 


III.2.1. Les méthodes de la théorie de la fonctionnelle de la densité (DFT) 


La théorie de la fonctionnelle de la densité (DFT) se distingue des méthodes ab initio basées 
sur les équations HF car elle s'appuie sur la notion de densité électronique p(r) et non sur celle de 
fonction d'onde multiélectronique. Cette approche est basée sur deux théorèmes dus à Hohenberg 
et Kohn [11]. 

Le premier stipule que "l'énergie totale d'un système dans son état fondamental ne dépend 
que de sa seule densité électronique p(r)" (l'énergie est donc une fonctionnelle de la densité 
électronique). Toute propriété de l’état fondamental, dont l’énergie, peut s’écrire comme une 


fonctionnelle de la densité électronique. 


E=E|p(r)| (IIL.1) 
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Le second théorème est l’analogue du principe variationnel appliqué à la densité et montre 


que la densité est stationnaire pour l’état fondamental. 
IIL.2.1.1. Les types de fonctionnelles existent 


e Les fonctionnelles locales (LDA) qui dépendent uniquement de la densité électronique en 
chaque point du système et négligent toute influence de l'inhomogénéité du système. 
Cette approximation est correcte lorsque la densité varie suffisamment lentement. Ces 
fonctionnelles (S [2], Xa [21], VWN [3]) ne sont plus utilisées à présent de façon 
régulière. 

e Les fonctionnelles à correction de gradient (GGA) qui introduisent dans leur expression le 
gradient de la densité permettent de tenir compte de l'inhomogénéité de la distribution 
électronique. Les fonctionnelles les plus employées sont : Becke88 [4], PWOI [5], P86 
[6] et LYP [7]. 

e Les fonctionnelles hybrides HF-DFT, apparues récemment, qui incluent pour l'énergie 
d'échange un mélange Hartree-Fock et DFT tandis que l'énergie de corrélation reste 


purement DFT. Ces fonctionnelles apparaissent comme étant les plus fiables du moment. 
a) La fonctionnelle B3LYP [8] 


B3LYP représente la fonctionnelle hybride la plus populaire de la chimie computationnelle. 
Elle a été introduite par l’équipe de Becke en 1993 [9]. La particularité de cette fonctionnelle est 
de présenter une combinaison linéaire entre des fonctionnelles d’échange-corrélation GGA et de 
l’échange Hartree-Fock. L’énergie d’échange calculée par cette méthode est composée de 80% 
de DFT et de 20% de HF. L’énergie totale d’échange corrélation peut être représentée par 


l’équation suivante : 


ER = EN + a (Et -E*)+a, (EE )+a (EC EPA) (UIL.2) 


Les indices x et c désignent l’énergie d’échange et de corrélation respectivement. LDA et 
GGA désignent les termes énergétiques calculés par la DFT. HF désigne la contribution calculée 
par la théorie HF. ao, ax et a. sont des coefficients constants définis empiriquement fixant le 
poids de chaque terme (0.2 ; 0.72 et 0.81 respectivement). 

Néanmoins, cette méthode présente quelques limitations telles que la sous-estimation des 
hauteurs de barrières énergétiques [10] et l’absence de prise en compte des interactions non- 
covalentes : la B3LYP est incapable de décrire des liaisons de Van Der Waals pour les composés 


liés par des interactions de portée moyenne. Malgré ces problèmes, cette fonctionnelle reste la 
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base de calculs pour la plupart des composés chimiques et l’outil le plus utilisé en modélisation 


moléculaire. 
b) La fonctionnelle B2PLYP 


Les fonctionnelles de densité double hybride (E?2”°?°) [11] sont considérées comme des 
méthodes de DFT «à rang supérieur» comparées aux fonctionnelles hybrides et méta-GGA. Ils 
mélangent un échange DFT standard (ÆE?"7) et une corrélation (E?°7) avec un échange HF 


(E#°) et une contribution supplémentaire à la théorie de la perturbation du second ordre ( E!7?) : 
ER" ={(1-a,)E," +a,Ef +(1-a )EPT a ET (IIL.3) 


La contribution ET? est obtenue par un terme de perturbation de Mgller-Plesset (MP2) basé sur 


les orbitales de Kohn-Sham optimisées de manière auto-cohérente par rapport aux trois premiers 


termes : 


gere = 1 5 5 [Gal ie) = Gb] ja) 


C 
4 ere; e, 6, 


(IL.4) 
ia jb 

L'exemple archétypal est la fonctionnelle B2-PLYP [12] qui utilise les fonctionnelles B88 
et LYP pour l'échange et la corrélation, en mélangeant l'échange HF (aX = 0,53) et la correction 
de perturbation (aC = 0,27). Le pourcentage d'échange HF dans les DHDF est généralement 


considérablement plus élevé que dans la plupart des fonctionnelles hybrides standards. Les 


DHDF 
E xC 


ont montré des performances robustes et avec des résultats généralement supérieurs à 
ceux d’autres méthodes de DFT dans de nombreux tests impliquant différents types 
d’énergétique [11,13-17]. Les études impliquant des systèmes à couplage d’échange sont 
relativement rares. Schwabe et Grimme [18] ont évalué la DFT double hybride sur un ensemble 
de composés organiques et quelques complexes de cuivre simples [19], confirmant ainsi la bonne 
performance de la méthode. Rajaraman et ses collaborateurs [20] ont rapporté que la 
fonctionnalité B2-PLYP à double hybride donnait de meilleurs résultats que les fonctionnelles 
hybrides tels que B3LYP, PBEO ou TPSSh pour la description du couplage magnétique dans les 
complexes Gd (I) -radicaux. Vogiatzis et al. [21] ont également signalé que B2-PLYP 
fonctionnait bien dans l'approche BS-DET pour les complexes à aubes dinucléaires, mais ils ont 


également noté qu'il était inférieur à B3LYP dans le cas d'un système à trichrome [21]. 
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c) La DFT-D3 

La méthode DFT-D3 présente en outre les très belles propriétés de complexité numérique 
mineure, que la dispersion intramoléculaire thermochimique «normale» est également améliorée 
de manière significative par rapport à de nombreuses fonctionnelles standard «aveugle de 
dispersion» et que les résultats peuvent facilement être analysés. Parce qu’elle combine les 
meilleures propriétés de toutes les autres approches, on peut s’attendre à ce qu’une sorte de DFT- 
D reste l’approche la plus largement utilisée du problème de la dispersion pendant au moins les 5 


à 10 prochaines années. 


IIL.2.1.2. Méthodes de la théorie de la fonctionnelle de la densité dépendante du temps (TD- 
DFT) 


La théorie de la fonctionnelle de la densité dépendante du temps (TD-DFT) est une approche 
qui décrit les excitations électroniques. Similairement au premier théorème de Hohenberg et 
Kohn [1] de la DFT, le théorème de Runge-Gross [22] établit une relation directe entre la densité 
électronique dépendante du temps et le potentiel dépendant du temps v /p(r, 4]; dès lors le 
potentiel et la fonction d’onde seront dépendants de la densité. 

Récemment van Leeuwen a généralisé le théorème de Runge-Gross aux systèmes 
moléculaires [23]. De la même manière, les équations de Kohn-Sham dépendantes du temps 
devront être dérivées. En supposant qu’il existe un système de référence dépendant du temps 
sans interaction avec un potentiel extérieur v,(r, t) dont la densité électronique p;(7, t) est égale à 


la densité exacte du système réel avec interaction p(F, t). 


N s) 
p(r,t) =p,(r,t) = De Gb] (UL.S) 


S1 le potentiel extérieur existe, on peut déduire les équations de Kohn et Sham dépendantes du 


temps : 
E[p@),t]= T[p(r)]+ [var [ptr, 0 [p(r, Ddr (UIL.6) 
pa',o .. 0Exc [P] 
veff [p@,t)|=v,.(t)+ dr'+ IIL.7 
CODAGE er TS (IIL.7) 
Il suffit ensuite de résoudre l’équation de schrôdiger dépendante du temps : 
. OP. (r,t 
[T, À Ver ] (rt) =i re (IIL.8) 
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Le calcul des états excités avec la TD-DFT est plus rapide que les méthodes ab initio. La 
théorie de la réponse linéaire est appliquée à la densité de charge d’une molécule soumise à un 
champ électrique oscillant dans le temps. La polarisabilité dynamique est obtenue à partir de la 
réponse du moment dipolaire, et les énergies d’excitations sont alors calculées comme les pôles 
et les forces d’oscillateurs comme les résidus de cette réponse, développés dans la base des états 
propres de l’opérateur hamiltonien DET. Des fonctions d’onde sont déterminées à posteriori pour 
les états excités, sous la forme de combinaisons linéaires de simples excitations à partir du 
déterminant Kohn-Sham. Cette méthode n’est vraie que dans la limite d’une fonctionnelle de la 
densité exacte. En pratique, elle est très sensible au mauvais comportement du potentiel 
d’échange-corrélation à longue portée, là où la densité est la plus polarisable. En conséquence les 
énergies calculées pour les états excités à transfert de charge peuvent s’en trouver gravement 
affectées [24, 25]. Cette méthode, adaptée aux systèmes de grande taille du fait de son coût 


modique, n’est donc pas encore fiable dans tous les cas, et doit être utilisée avec précaution. 


1I1.2.1.3. Solvatation 


La majorité des réactions chimiques et biologiques ont lieu en solution et les effets dus au 
solvant peuvent être très importants. Celui-ci influe directement sur les interactions non liantes et 
sa présence peut être simulée de deux manières : en introduisant effectivement les molécules de 
solvant dans le système étudié (solvant explicite) ou bien en mimant sa présence implicitement 


dans les calculs de la fonction d’énergie (solvant implicite). 


IIL.2.1.4. L’analyse Orbitalaire NBO 


L’idée originale de la méthode NBO est de proposer un découpage de la densité électronique 
d’un système moléculaire pour décrire la structure électronique de la molécule dans un schéma 
de type Lewis [26]. Ainsi la densité est décrite sur et entre les atomes avec des doublets de cœur 
(CR), des doublets non-liants (LP), des doublets de liaison (BD) résultants de la mise en 
commun par deux atomes de deux électrons dans des orbitales hybrides. En effet, la 
diagonalisation de la matrice densité d’ordre un conduit à l’obtention des orbitales naturelles qui 
sont les fonctions propres de cette matrice. Les valeurs propres représentent alors les occupations 
fractionnaires des orbitales naturelles. Cette approche permet de calculer correctement les 
diverses propriétés monoélectroniques. 

Inspiré par cette approche, Weinhold [26] a proposé l’analyse NBO dont la procédure se 
déroule en trois étapes. La première étape consiste à diagonaliser la matrice densité d’ordre un en 
blocs monocentriques atomiques sans diagonaliser la matrice dans son entier (orbitales 


atomiques naturelles NAO). La seconde étape du processus NBO consiste à orthogonaliser les 
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orbitales NAO. La dernière étape est une recherche d’orbitaux naturels hybrides, processus qui 
diagonalise des blocs bicentriques ou tricentriques d’orbitales naturelles orthogonalisées. Cette 
analyse est cohérente en général avec les structures de Lewis, car la troisième étape commence 
par sélectionner les blocs monocentriques avec une occupation >1, 90, ce qui correspond aux 
orbitales monocentriques de cœur ou de valences (paires libres et liaisons). Les autres orbitales 
sont ensuite orthogonalisées par blocs bicentriques. 

Dans l’espace NBO, chaque orbitale naturelle de liaison (doublement occupée : la structure 


naturelle de Lewis) est associée à une orbitale “antibonding” (non-Lewis) : 


Op = Cala +Cehg (IIL.9) 


Os = Cha +Cuhz (HL.10) 


Où hi est une orbitale naturelle hybride et ci son coefficient de polarisation de l’atome i. 

Une analyse perturbative des interactions “donneur-accepteur”, appelée “analyse perturbative 
au second-ordre” est effectuée entre les NBO occupées et vacantes. Elle mesure le transfert 
d’énergie engendré par la combinaison d’une NBO occupée (i d’énergie E;) et d’une NBO 
vacante (j d’énergie E;) en prenant en compte l’élément hors-diagonal de la matrice de Fock F (i ; 


j) ainsi que l’occupation qj; de la NBO :: 


E,=AB;=Q ee Fu CD (IL.11) 


F? ul 
Les avantages de la théorie NBO 
Cette techniques de localisation est très utilisée dans la chimie théorique moderne, elle permet 

de valider les structures de Lewis pour des systèmes avec “ Bonne”” liaisons et elle peut être 
aussi serviable dans l’analyse des systèmes avec des liaisons délocalisées. L’analyse présente les 
propriétés suivante : 

Ÿ”_ Applicable à tout niveau de calcul. 

Ÿ__ Un bon fonctionnement computationnel. 

Ÿ”_ Construit la meilleure structure de Lewis. 


Ÿ”_ Applicable pour une variété de liaisons. 


IIL.2.1.5. Les descripteurs de réactivité globale et les indices de Fukui 


Les paramètres de la réactivité globale sont utilisés pour reconnaître le lien entre la structure, 
la stabilité et la réactivité chimique globale de toutes les molécules. Ces paramètres sont utilisés 
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dans la progression des relations quantitatives structure — propriété, structure — activité et 
structure — toxicité. Et il est bien connu qu’il existe une relation entre aromaticité et dureté [27]. 
Dans les recherches actuelles, nous avons utilisé la DFT qui fournit des éclaircissements sur les 
perceptions universelles significatives de la stabilité et de la réactivité de la structure moléculaire 
[28]. Actuellement, les descripteurs de réactivité chimique globaux tels que la dureté (n), le 
potentiel chimique (li), la douceur (S), l'électronégativité (x) et l'indice d'électrophilicité (w) du 
composé du titre ont été calculés à l'aide des valeurs EHOMO et ELUMO sous forme 
d'ionisation potentiel (1) et affinité électronique (1) respectivement. Ces paramètres du composé 
indiqué ont été calculés comme suit : En utilisant le théorème de Koopman [291 pour les 
molécules à enveloppe fermée, la dureté chimique de toute molécule peut être calculée par la 


relation suivante. 
n = (IL.12) 


Le potentiel chimique calculé par : 


u={04) (11.13) 


Le softness calculé par : 


S=— II.14 
2 QI. 14) 


L’électronégativité calculé par : 


r-[54) (IL.15) 


L’électrophilicité calculé par : 
nus (HIL.16) 


Des paramètres de réactivité locaux sont nécessaires pour différencier différents sites réactifs 
d'une molécule. Parr et Yang ont montré que chez les espèces chimiques, les sites avec les 
valeurs les plus élevées de la fonction de Fukui (f (r)) sont ceux avec une réactivité plus élevée. 


La fonction de Fukui est définie comme [30] : 
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f= 2 | (HL.17) 
ON |, 

Où f (r) est la densité électronique, N le nombre d'électrons et v le potentiel externe exercé par 
le noyau. La fonction de Fukui est un descripteur de réactivité local qui indique les régions 
préférées où une espèce chimique va changer de densité lorsque le nombre d'électrons est 
modifié. Par conséquent, il indique la propension de la densité électronique à se déformer à une 
position donnée lors de l'acceptation ou du don d'électrons [31-32]. Pour étudier la réactivité au 
niveau atomique, il convient d'utiliser un moyen plus pratique de calculer les fonctions f (r) à 
résolution atomique, les fonctions de Fukui à atome condensé. Pour un atome K dans une 


molécule : 


fi = qu N+D-q(N) 
fe = qi N)-q,(N -D (IIL.18) 


0 1 
d = LAN +D—a (ND) 


Où fk *, fc et fi° sont respectivement des attaques nucléophiles, électrophiles ou radicales et qk 
est la population électronique de l'atome Kk de la molécule considérée évaluée à partir de l'analyse 
de population de Mulliken. La plus grande valeur de la fonction Fukui implique la plus grande 


réactivité du site correspondant [33]. 


IIL.2.1.6. La surface du potentiel électrostatique moléculaire (MEP) 


La surface du potentiel électrostatique moléculaire (MEP) est une méthode de cartographie du 
potentiel électrostatique sur la surface de densité isoélectronique. MEP affiche simultanément la 
taille moléculaire, la forme et les régions de potentiel électrostatique en termes de gradation de 
couleur et constitue un outil très utile dans la recherche de corrélation entre la structure 
moléculaire et la relation de propriété physicochimique de molécules comprenant des 
biomolécules et des médicaments [34-39]. La palette de couleurs de la surface du MEP est la 
suivante : rouge pour une charge riche en électrons, partiellement négative ; bleu pour les 
charges partiellement positives déficientes en électrons ; bleu clair pour la région légèrement 
déficiente en électrons ; jaune pour les régions légèrement riches en électrons ; vert pour neutre 


(potentiel zéro) ; respectivement [40]. 
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IIL.2.1.7. Propriétés optiques non linéaires (NLO) 


Les matériaux à forte réponse d'optique non linéaire (NLO) sont très utiles dans les dispositifs 
optoélectroniques, et l'optique non linéaire trouve de grandes applications dans les technologies 
de l'information et dans d'autres industries [41]. Les matériaux NLO ont également une large 
gamme d'applications pour les dispositifs de communication photonique et de mémoire 
numérique, la défense nationale et l'industrie pharmaceutique [42]. Les composés ayant des 
groupes donneurs d'électrons et accepteurs d'électrons avec un système à conjugaison x peuvent 
être considérés comme des candidats puissants pour de telles applications, car l'augmentation de 
la conjugaison dans la molécule entraîne une modification des propriétés NLO lorsque les 
propriétés NLO d'une molécule sont liées au fossé énergétique entre HOMO et LUMO. La 
valeur des paramètres ci-dessus de tous les composés a été simulée au niveau théorique de 
B3LYP / 6-311++G (d, p) avec les calculs supplémentaires calculés en utilisant les équations 
suivantes : 


Le moment dipolaire a été calculé par : 


1/2 


Un = + Hs HE) (L.19) 
La polarisabilité totale a été calculée par : 

1 
Lio 7 (a, +a,, +@..) (HI.20) 
L'anisotropie de la polarisabilité a été calculée par : 
AG = F [a -a,) +(a, -a,) +(a, -a,.) +642] (L21) 
Les composants de la première hyper polarisabilité peuvent être calculés par : 
B=p: Ds (UIL.22) 


En utilisant les composantes x, y, z, la magnitude de la première hyperpolarisabilité (ft) peut 


être calculée comme suit : 


Ba = CB? + B2 + 87) (IIL.23) 


Où PB, By et B, sont Bx = (Brx + Brxy + Bryy), BY = (Byyy + Baxz + Byyz), BZ = (Bxez + Byzz + Pare). 
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Par conséquent, l'équation finale pour la magnitude du premier calcul total de 


l’hyperpolarisabilité statique totale est donnée par : 


1 
Pi = [Bsxx + Pi + Be + (ee LA BE pu ps) + (Be + Bas + B,) 1 (UL.24) 
IIL.2.1.8. Choix des méthodes de calcul 


Les calculs de la modélisation moléculaire ont été effectués au moyen du logiciel 
GAUSSIAN 09 [43]. Avec le même ensemble de bases 6-311 ++ g (d, p), les calculs ont été 
effectués au niveau DFT en utilisant quatre fonctionnelles d’échange : la fonctionnelle hybride 
(méthode hybride à trois paramètres de Becke et fonctionnelle B3LYP de corrélation Lee — Yang 
— Parr) [44]), le double hybride (B2PLYP [45]) et leurs versions de dispersion Grimme : 
B3LYP-D3 (BJ) et B2PLYP-D3 respectivement. 

Dans le modèle de calcul des complexes, l'anion était ignoré et le complexe cationique était 
pris en compte. Le calcul de l'état excité par TD-DFT [46,47] a été effectué au même niveau 
théorique, le méthanol a été utilisé comme solvant et le modèle de type polarisable en continu 
CPCM [48] semblable à un conducteur a été appliqué comme modèle de solvatation. 

Les structures optimisées ont été vérifiées comme des minima sur les surfaces d'énergie 
potentielle par des calculs de fréquence. Les orbitales des frontières moléculaires HOMO 
(orbitale moléculaire la plus occupée), LUMO (orbitale moléculaire la plus basse inoccupée) et 


les structures optimisées ont été visualisées avec GaussView [491]. 


IIL.3. Docking moléculaire 


Le docking moléculaire consiste a déterminé la position relative de deux molécules (un ligand 
et un récepteur). La structure obtenue confère les propriétés à l’ensemble (le complexe) ainsi 
formé. Par exemple, la recherche de bonne structure est cruciale dans la conception de nouveau 
médicament. Ou encore, selon l’association de 2 protéines, le signal déclenché lors de leur 
association peut être différent. D’un point de vue informatique, la prédiction de structure se 
traduit par un problème d’optimisation. À chaque position relative possible est associé une 
énergie (également appelée score, fitness, etc.). La structure la plus probable est alors celle qui 


minimise cette énergie. 
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Complexe + AG 
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Docking Scoring 


Figure 3.1 : Principe générale d’un programme de docking. 


Les interactions protéine-protéine et protéine-ligand jouent un rôle clé dans 
l’organisation des systèmes biologiques. Elles permettent la régulation de certains processus 
biologiques [50] ainsi que la transmission des signaux ou encore la catalyse de diverses 
réactions biochimiques. Connaître la façon dont les protéines interagissent avec d’autres 
entités biochimiques est une étape essentielle en vue de comprendre les processus biologiques 
dans lesquels elles sont impliquées. Le développement d’approches prédictives ouvre la voie 
à la conception assistée par ordinateur de systèmes protéiniques aux propriétés modifiées et 
présente donc un intérêt indéniable pour la recherche et l’industrie pharmaceutique et 
médicale. 

Dans le domaine de la chimie théorique, le docking moléculaire est une méthode qui 
prédit la conformation (position et orientation relatives) la plus favorable de deux molécules 
en interaction et formant un complexe stable. La connaissance de cette conformation 
préférentielle permet par la suite l’estimation de la force d’association (affinité de liaison) 
entre ces deux molécules. Cependant, chaque programme de docking diffère par rapport aux 
autres par son algorithme d’échantillonnage, sa manière de manipuler la flexibilité du ligand 
et de la protéine ainsi que sa fonction d’évaluation des complexes. Il existe deux grands types 
de docking moléculaire. Le docking rigide consiste à obtenir la conformation préférentielle 
d’un système protéine-ligand en considérant que chacune des deux molécules conserve une 
géométrie interne fixe. Dans ce cas, la relaxation de la géométrie interne de chaque entité, en 
interaction dans le complexe, n’est pas prise en compte [50]. Actuellement cette approche est 
utilisée dans le docking protéine-protéine. L’approche la plus récente dans le traitement de la 
flexibilité, considère la cible protéique comme corps rigide tandis que le ligand est flexible. 
En effet elle représente l’approche la plus utilisée dans la plupart des logiciels ou programmes 


de docking. Néanmoins, celle-ci reste approximative, du fait que le récepteur ne doit pas être 
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considéré comme corps rigide (au moins la flexibilité du site d’interaction doit être prise en 


compte). 
IIL.3.1. Principe général Docking moléculaire protéine/ligand 


Le docking (amarrage) moléculaire est une méthode empirique qui permet de prédire 
l’affinité entre deux molécules. Il est généralement utilisé pour prédire l’affinité d’un ligand pour 
une protéine, entre autre, la position et l’orientation la plus favorable pour un ligand venant en 
interaction avec une protéine cible. Cet outil est actuellement très sollicité dans la recherche de 
nouvelles molécules à visée thérapeutique. Il est utilisé pour faire des criblages très rapides et 
peu coûteux afin d’identifier une molécule pouvant avoir une forte activité sur une protéine cible 
ou au contraire d’identifier la cible d’une molécule ayant des activités biologiques intéressantes. 
D’autres parts, lorsque la structure tridimensionnelle d’une protéine cible est définie, le docking 
permet de mettre en évidence les régions et les résidus à explorer pour optimiser l’affinité d’un 
ligand avec cette cible. Aïnsi, 1l est également utilisé pour optimiser la sélectivité d’une molécule 
entre deux ou plusieurs protéines. Le principe général est, considérant une cible de géométrie 
donnée, d’essayer de placer le ligand en interaction avec la protéine de la façon la plus favorable 
possible. Lors du docking, la protéine cible est généralement figée alors qu’une certaine 
flexibilité est attribuée au ligand (angle de torsion). Il existe de nombreux logiciels de docking, 
qu’on peut répartir en deux catégories selon la méthode utilisée : la simulation de trajectoire et 
l’ajustement de fragments. A partir d’une position initiale aléatoire, la simulation de trajectoire 
consiste à explorer plusieurs positions (avec différentes conformations du ligand) par évaluation 
de l’énergie d’interaction ligand/protéine. La meilleure position est celle qui est la plus basse en 
énergie, donc correspond au complexe ligand/protéine le plus stable. La seconde approche, 
ajustement de fragments, est basée sur la complémentarité des formes. Dans cette méthode, le 
ligand est décomposé en plusieurs fragments puis reconstruit dans le site actif, en essayant de 
faire correspondre les géométries et les fonctions chimiques. Cette méthode est beaucoup plus 
rapide que la simulation de trajectoire et permet généralement le criblage rapide de vastes 
librairies de molécules. La méthode par simulation de trajectoire est plus lente mais prend mieux 
en compte la flexibilité du ligand et permet l’exploration de régions plus vastes. Compte tenu des 
avantages et inconvémients de chaque méthode, deux logiciels seront utilisés dans cette thèse : 
AUTODOCK 3.05 (la simulation de trajectoire), pour la validation du modèle construit par 


homologie de séquence et GOLD (ajustement de fragments), pour les études de prédiction. 
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IIL.3.2. Docking avec AUTODOCK 


L’interaction ligand/protéine est évaluée en suivant plusieurs étapes. 


IIL.3.2.1. Préparation des sources : 


Il s’agit de vérifier et d’affiner la structure du ligand et de la cible (protéine) avant de les 
sauvegarder au bon format mol2. Tous les atomes d’hydrogène du ligand sont présents, alors que 
seuls les hydrogènes essentiels (polaires) de la cible ne sont gardés. Les charges partielles sont 
assignées par la méthode qui convient (Gasteiger pour le ligand et KOLL UNI pour la cible). 
Cette étape est effectuée avec le logiciel SYBYL et son module BIDPOLYMER. Le fichier 
cible.mol2 est converti en format pdbq puis pdbgt avec le module convpdbq d’AUTODOCK 
dans lesquels on trouve les charges et paramètres de solvatation des atomes de la cible. Pour le 
ligand.mol?, on définit avec le module deftors d'AUTODOCK les pivots qui sont en libre 


rotation durant la simulation et on génère le fichier ligand.pdbq. 
IIL.3.2.2. Génération d’une boîte contenant la grille et calcul des grilles de potentiels : 


Les deux fichiers ligand.pdbqt et cible.pdpqt ayant des coordonnées spatiales différentes, on 
rapproche le ligand de la cible (site catalytique). Une boîte contenant une grille est définie autour 
de la cible (cible.pdbqf), qui correspond à la région où peut s’amarrer le ligand (/igand.pdbgf). 
Cette boîte est généralement centrée sur le ligand et ses dimensions doivent être éditées afin de 
l’ajuster au mieux à la région de la cible qu’on souhaite étudier (site catalytique). Elle est 
générée par AUTODOCK et visualisée sur SYBYL (Figure 3.2). Pour chaque type d’atome 
présent dans le ligand (C, O, N...), le programme AUTOGRID de AUTODOCK calcule 
l’énergie d’interaction avec la cible lorsqu'il est placé en un point donné de la grille. Cette 
évaluation est répétée pour chaque type d’atome et chaque point de la grille. Le but est d’évaluer 


très rapidement l’énergie pour les différentes configurations du ligand. 
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Figure 3.2 : Une boîte contenant la grille par docking moléculaire avec AUTODOCK 
IIL.3.2.3. Génération des paramètres de recherche : 


On génère le fichier ligand.cible.dpf à partir du ligand.pdbg et cible.pdbgs. Ce fichier qui 
contient les paramètres de recherche est par la suite édité pour ajuster les paramètres tels que le 


nombre de conformations générés, la précision de recherche et bien d’autres. 


IIL.3.2.4. Recherche et propositions de solutions par AUTODOCK : 


Le programme procède à la recherche du complexe ligand/cible le plus favorable en déplaçant 
et en changeant la conformation du ligand. La recherche est faite via un algorithme génétique qui 
permet de générer des configurations et positions du ligand rapidement. À chaque configuration, 
l’énergie d’interaction est évaluée en mécanique moléculaire avec un champ de force spécifique 


à AUTODOCK et en utilisant la grille calculée précédemment. 


IIL.3.2.5. Analyse de résultats : 


A la fin du docking, un fichier ligand.cible.dlg est généré. Il contient des histogrammes (ou 
clusters) regroupant les solutions proposées par similitude et les classe par énergies croissantes 
(Figure 3.3). La similarité entre deux solutions est calculée par la racine de la moyenne des 
carrées des écarts (Root Mean Square Deviation-RMSD). Deux solutions classées dans le même 
cluster, si le RMSD entre les deux est inférieur à une valeur seuil. Cette valeur est définie lors de 


la génération des paramètres, dans le fichier ligand.cible.dpf. Théoriquement, la meilleure 
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solution est celle qui a l’énergie la plus basse et se trouve idéalement dans le cluster le plus 
peuplé. 
On peut ensuite extraire les coordonnées des solutions trouvées pour les visualiser avec des 


logiciels de visualisation (VMD, PYMOL, SYBYL et d’autres). 
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Figure 3.3 : Résultat de docking : Histogramme par Docking moléculaire avec AUTODOCK 


IIL.3.3. Docking avec GOLD : 


GOLD (Genetic Optimisation for Ligand Docking) est un logiciel de docking basé sur la 
complémentarité de forme et il est généralement utilisé pour le docking de petites molécules 
organiques dans le site catalytique d’une protéine. En comparaison à AUTODOCK, GOLD est 
plus rapide et très simple d’utilisation. Le logiciel est fourni sous forme d’un ensemble de 
programme (HERMES) pour la visualisation et la manipulation des structures. Pour faire un 
docking, on ouvre HERMES et on rentre dans le programme GOLD. La première étape consiste 
à charger la protéine cible sous format mol2 et de vérifier que tous les atomes d’hydrogène sont 
présents. On définit ensuite une région à étudier autour de la cible, qui correspond au site 
catalytique de la protéine. Dans la deuxième étape, comme pour la cible, on charge le ligand sous 
format mol2 et on vérifie que tous les atomes d’hydrogène sont présents. A cette étape, on définit 
le nombre de cycle (GA runs) qui correspond au nombre de solutions générées à la fin du 
docking. Contrairement à la cible qui est figée, on attribue une certaine flexibilité au ligand. La 
troisième étape consiste à définir une fonction de score, qui permet d’évaluer l’affinité entre le 
ligand et la cible. Il existe plusieurs fonctions de score et CHEMPLP, qui est la fonction de score 
par default de GOLD est celle qui a été utilisée dans cette thèse. C’est une fonction de score 
empirique qui tient compte de la complémentarité stérique entre le ligand et la cible, et inclut des 
termes de distance et d’angle pour les liaisons hydrogènes. Dans la dernière étape, GOLD 


procède à la recherche de la position du ligand qui exploite le mieux le site catalytique de la cible 
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via un algorithme génétique. A la fin du docking, plusieurs solutions sont générées en fonction 
du nombre de cycle (GA runs) défini. Ces solutions sont classées en fonction de leur score et la 
meilleure pose de docking correspond à celle qui a le score le plus haut (best ranking) et qui 


revient le plus de fois. 


IIL.3.4. Le scoring 


La procédure de docking permet de générer une liste de complexes représentant les 
modes d’association favorables entre le ligand et le récepteur macromoléculaire. Les 
méthodes de "scoring" permettent, par la suite, d’évaluer l’énergie de liaison du complexe 
formé et de donner un score aux poses obtenues lors de la phase de docking, afin de trouver 
celui ou ceux les plus susceptibles de reproduire aux mieux le mode d’association réel [51]. 
L'association non covalente et réversible d'un récepteur (R) et d'un ligand (L) pour former un 
complexe récepteur-ligand (RL), dans un milieu aqueux contenant des électrolytes peut être 


décrite par l'équation (25) : 
R+LSRL (IL.25) 


À l'équilibre, cette réaction est déterminée par un constant d’équilibre Key (équation 26) 


et aussi par l'énergie libre de liaison AG;,,, (équation 27). 
K,,= [RL] /[R]x[£] (UIL26) 
AGjna = — RTINK,, (IL27) 


Où R est la constante des gaz parfaits et T est la température (en kelvin) [52, 53]. 
Théoriquement, le complexe est favorable si la variation d’énergie libre globale de 
complexation est négative (AGbind de complexation < 0). En pratique, l’évaluation de l’énergie 
libre des complexes est une tache souvent couteuse ce qui limite son utilisation en routine. De 
ce fait, des méthodes approximatives ont été développées pour distinguer (évaluer et classer) 
les meilleurs complexes parmi ceux générés par une procédure de docking : les fonctions de 
score. Dans les programmes de docking, on trouve différents types de fonction de score : 
celles utilisant un champ de force de mécanique moléculaire (CHARMM, AMBER, GAFF, 
GROMOS, OPLS, MMFF..etc.), celles reposant sur les connaissances actuelles 
(Knowledge-based) et les méthodes empiriques (Bohm's Score, ChemScore, Fresno, F-Score, 
SCORE, X-Score....etc.). Les fonctions de score les plus couramment utilisées sont présentées 


schématiquement dans la Figure 3.4 [$1]. Suivant la structure du ligand et son emplacement vis- 
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à-vis du récepteur, AutoDockné value leur affinité en utilisant une fonction de score basée sur 


une forme empirique de l'énergie libre d'interaction entre le ligand et le récepteur [51-54]. 


-AutoDock Score Fonction 


Utilisant un champ de -GoldScore 
+ | force de mécanique ———+ -G-Score 
moléculiare -D-Score 
——— -DrugScore 
Knowledge-based + -ITScore 
s -PMF 
-LigScore 
SCORE, X-Score 
Empiriques +  F-Score et ChemScore 


Figure 3.4 : Différents types de fonctions de score les plus connues (en vert) et les 


exemples (en bleu) les plus fréquemment cités dans la littérature [511]. 


L'énergie est donnée par l’équation (28) : 


AG nd 5 AG + AG, pond 7 AG, 


elec 


+ AG. + AG., (1.28) 


Les trois premiers termes sont des termes classiques de mécanique moléculaire ; il s’agit 
respectivement des énergies : dispersion/répulsion des atomes, liaisons hydrogène et 
interactions électrostatiques. AGtr est un terme qui traduit l’augmentation d’énergie du 
système due à la restriction des rotors libres du ligand et à la restriction des rotations et 
translations du ligand lors de la complexation au récepteur. Cette perte de degrés de liberté est 
une perte d’entropie. AGsai est un autre terme lié à l’entropie qui décrit les variations d’énergie 
du système lors de la désolvatation du ligand au moment de la complexation au récepteur. 
C’est une modélisation partielle de ce que l’on appelle l’effet hydrophobe dû aux variations 
de l’entropie du solvant aux interfaces solvant-soluté [55]. 


IIL.3.5. Prédiction des interactions par docking moléculaire 


L'utilisation des méthodes de docking (ancrage, amarrage ou arrimage moléculaire) 
dans le processus de conception de médicaments a débuté il y a plus de 30 ans [56]. Leur objectif 
était de prévoir la capacité ou non d’une molécule à se lier au site actif d’une protéine en se 


basant sur la prédiction de la conformation et de l’orientation de la molécule lors de sa liaison au 
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récepteur [57]. A cet effet, les méthodes de docking combinent l’utilisation d’un algorithme de 
recherche, permettant de générer des modes de liaisons putatives du ligand dans le récepteur, ou 
« poses », et d’une fonction de score, employée pour classer les différentes poses selon un score 
prédit d’affinité [58]. Les méthodes de docking s’attachent donc, d’une part, à identifier les 
molécules qui sont des ligands véritables du récepteur parmi toutes celles étudiées, mais aussi 
d’autre part, à déterminer les poses correctes soit les conformations adoptées par les ligands lors 


de la liaison au récepteur. 
1II.3.6. Les interactions intermoléculaires Protéine-Inhibiteurs 


Au niveau du site actif, le processus d’interaction se fait par liaison complémentaire du ligand 


avec l’architecture protéique. Cette association est assurée grâce à plusieurs types de liaisons 


faibles. 
11L1.3.6.1. Les différentes interactions 


e Les interactions électrostatiques : Elles agissent à plus grande distance entre dipôles chargés. 
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Figure 3.5 : Interactions électrostatiques 


e Les forces de Van Der Walls: Elles résultent de l’interaction des nuages électroniques de 
deux atomes adjacents conduisant à la présence d’une force attractive pour des distances de 
3-4 À. Les forces de van der Waals sont généralement faibles. Leur nombre élevé va 


produire au total une force importante [59,60]. 
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Figure 3.6 : Les forces de Van der Waals 
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e La liaison hydrogène : Cette liaison intervient lorsqu’un atome d’hydrogène lié à un atome 
électronégatif (le donneur) est attiré par un autre atome électronégatif (l’accepteur) et elle agit à 
très courte distance (0,8 à 2,8 À). Les liaisons hydrogènes sont peu nombreuses et s’adaptent très 
bien à la flexibilité (l’angle peut varier de 120° à 180°). Les liaisons hydrogènes sont définie par ; 
la distance entre les deux atomes qui forment la liaison covalente D-H de type à, la longueur de 


l’interaction H...A et l’angle D-H...A[59,601. 
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Figure 3.7 : Liaison hydrogène 


e Les interactions hydrophobes : Les molécules dépourvues de groupes chargés ou d'atomes 
capables de former des liaisons hydrogène ne peuvent donc pas s’hydrater ; 
pour cette raison on les nomme substances hydrophobes. L’effet hydrophobe est la tendance 


de ces groupes à se rassembler par coalescence de façon à minimiser les contacts avec l’eau. 


Figure 3.8 : Interactions hydrophobes 
IIL.3.6.2. Implication des interactions dans la structure et la flexibilité 


Les conditions environnementales (température, pH, salinité...) affectent la stabilité et la 
flexibilité des biomolécules. Leur rigidité structurale est suffisante pour préserver leur forme 
native spécifique ; leur flexibilité interne permet leur bon fonctionnement (activité 


enzymatique des enzymes...). La flexibilité est, donc, essentielle pour conserver leur activité 


biologique. 
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Les liaisons de faible énergie sont nécessaires à l’établissement et au maintien des 
conformations De même que pour le fonctionnement des mécanismes des processus 
biologiques. Leur grand nombre, par effet cumulatif, leur confère un rôle capital au niveau du 
repliement tridimensionnel et de la stabilité des structures. En revanche, la facilité quand ces 
liaisons faibles à se rompre permet les interactions moléculaires caractéristiques des réactions 
biochimiques dont découlent les propriétés de la vie : en permanence, les molécules se 
reconnaissent, s’assemblent, réagissent et se dissocient. D’autre part, elles permettent une grande 
souplesse conformationnelle, clé de la flexibilité des molécules biologiques, et donc des 


capacités d’adaptation, de reconnaissance et de spécificité [61]. 
IIL.3.6.3. Importance de la structure co-cristallisée pour un criblage virtuel 


Lorsque la structure 3D de la cible biologique d'intérêt est disponible, des méthodes 
dites « basées sur la structure » peuvent être employées pour réaliser le criblage virtuel. Il 
existe deux types de méthodes expérimentales pour obtenir la structure 3D d’une cible : la 
diffraction des rayons X et la résonnance magnétique nucléaire (RMN). Les méthodes de 
prédiction, telles que les méthodes de docking qui sont les plus populaires, utilisent la structure 
3D de la cible pour découvrir de nouveaux composés actifs. Ces techniques nécessitent 


l’identification préalable du site de liaison. 
IIL.3.6.3.1. Identification du site actif 


L'identification du site de liaison (ou site actif) est nécessaire et cruciale pour pouvoir 
réaliser un criblage virtuel basé sur la structure. Pour cela, une solution est possible : obtenir 


une structure 3D de la protéine co-cristallisée avec un ligand [621]. 
IIL.3.6.3.2. Structure co-cristallisée avec un ligand 


L’analyse de la structure cristallisée d’une protéine avec un ligand permet de rechercher 
les interactions s’établissant entre ces deux partenaires [63]. Les acides aminés clés, c’est-à-dire 
impliqués dans les interactions ligand-protéine, permettent ainsi de guider la définition du site 
de liaison. Il s’agit de la méthode la plus sûre pour identifier un site de liaison puisqu’elle 


utilise des données expérimentales et non pas des données prédites. 
IIL.3.6.4. Réalisation du criblage virtuel (Docking moléculaire) 


Les calculs de Docking moléculaire ont été essentiellement effectués pour prédire et 
comprendre les modes d’interactions de nos inhibiteurs, non disponibles dans la PDB (Protein 
Data Bank), à la cavité du site actif de deux isoformes cytosoliques, Echer et candidas A. 
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Dans cette étude, nous avons utilisé le programme AutoDock 4.2.6 dans sa version la plus 
récente. AutoDock est un ensemble de trois sous-programmes destinés pour prédire la 
conformation la plus favorable d’un ligand au sein de son récepteur. Le contrôle des paramètres 
définissant les propriétés du ligand est réalisé avec le sous-programme AutoTors. Le sous- 
programme AutoGrid effectue, de son côté, le calcul des différentes grilles de potentiel. Enfin, le 
sous-programme AutoDock réalise les calculs d'amarrage moléculaire proprement dits. Nous 
avons également utilisé l'interface graphique "ADT" (AutoDock Tools) version 1.5.6 pour 
visualiser la localisation précise de l'espace à l’intérieur duquel vont se faire les calculs [64]. La 


procédure de simulation par le programme AutoDock 4.2.6 se décompose en 5 étapes : 
IIL.3.6.4.1. Préparation des molécules à l’arrimage 


Le complexe protéine-Inhibiteur est téléchargé dans le format pdb à partir de la banque 
de données en introduisant son code ID. Parmi les nombreux complexes proposés par la PDB, 
nous avons choisi les codes 1H1W et 1HNJ, car ils s’agissent des structures correctement 
définies avec des résolutions de 1.55 À et 1.98 À Avant leur utilisation par AutoDock, les 
complexes 1H1WI[65] et 1HNJ[66] sont séparés de leurs ligands afin d’obtenir des récepteurs 
avec des sites actifs libres. Chaque récepteur est ensuite préparé en définissant les charges 
partielles de ses atomes par l'intermédiaire du sous-programme AutoTors qui attribue des 
charges de type "Kollman” calculées avec le champ de forces empirique Amber [67].IL est aussi 
nécessaire de rajouter, aux récepteurs, les hydrogènes dits polaires c’est à dire « essentiels ». Les 
fichiers des protéines ainsi préparées sont convertis par AutoDock au format pdbq qui 
correspond aux données des charges. Enfin, on intègre dans ces fichiers les paramètres de 
solvatation des atomes avant de l’enregistrer dans le format pdbqt (format propre à 
AutoDock).Les géométries initiales des ligands utilisés dans ce travail sont optimisées 
précédemment (voir le paragraphe II du chapitre Il) avec la méthode DEFT pour minimiser leur 
énergie. Les molécules ainsi obtenues sont enregistrées au format pdb. La phase finale de 
préparation du ligand est réalisée par le sous-programme AutoTors4.Ce dernier a pour rôle de 
contrôler les mouvements du ligand, en définissant les atomes aromatiques, les liaisons flexibles 
et les angles de torsion. La molécule finale est enregistrée dans le format pdbqt avant son 


utilisation par AutoDock 4.2.6. 
IIL.3.6.4.2. Calcul des grilles de potentiel 


Pour augmenter la vitesse d’évaluation de l’énergie du système, le récepteur (soit ACHI ou Il) 
est plongé dans une grille tridimensionnelle englobant largement le site actif de la 
protéine étudiée et permettant la libre rotation du ligand dans ce site. Le centre de cette boite est 
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déterminé par les coordonnées X, Y et Z de -42.806, 19.079 et 2.962 pour l'isoforme Escherichia 
coli et de -25,208, 16,643 et 5,087 dans le cas de candidas, respectivement, avec les dimensions 
50*50*50 À3 (Figure 3.9). L’espacement ou la maille de la grille est fixé par défaut à 0.375 À 
(correspondant au quart de la longueur d’une simple liaison C-C). La boîte est alors centrée sur 
le ligand dont les dimensions sont proportionnelles à la taille de tous les ligands étudiés. L’étape 
suivante du programme AutoDock 4.2.6 consiste à exécuter les instructions grâce au sous- 


programme AutoGrid4. Il s’agit de calculer l’affinité entre les différents types d’atomes, de la 


protéine (C, H, O, N et S) et le ligand, à l'intérieur de la boîte. 


Figure 3.9 : Position de la boîte d’amarrage dans le site actif de 1H1W 
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Chapitre IV 


Résultats et discussions 


IV.1. Introduction 


Le présent chapitre qui est réservé principalement aux résultats et discussion, comporte au 
début toutes les descriptions et les conditions du travail. En effet, après une présentation des 
techniques et des appareillages utilisés, nous décrirons les modes opératoires adoptés pour 
l’élaboration des ligands ainsi que de leurs complexes de cobalt (ID et cuivre (I). Ces complexes 
ont été étudiés en phases liquide et solide. Les études en solution ont été réalisées par 
spectrophotométrie UV-vis à 25°C et les résultats obtenus ont été traités pour obtenir les 
constantes de stabilité et les distributions des espèces. A l’état solide, les complexes synthétisés 
ont été caractérisés par différentes méthodes telle que FT-IR, RMN 'H et par ATG/DSC. 

Aujourd’hui, la recherche et la synthèse de nouveaux composés chimiques sont, souvent, 
associées à une étude par modélisation moléculaire [1]qui constitue une approche scientifique 
complémentaire permettant de travailler sur des sujets extrêmement variés, que ce soit dans un 
but fondamental ou beaucoup plus appliqué. Dans ce travail, nous avons utilisé la modélisation 
moléculaire par des méthodes basées sur la DFT (Density Functional Theory). 

Le principe d’un calcul consiste à minimiser l’énergie d’une molécule, décrite dans le vide ou 
dans un solvant, à partir de la position dans l'espace des atomes qui la constituent. La structure 
tridimensionnelle de plus basse énergie, appelée modèle, peut être considérée comme une 
représentation de la conformation de la molécule "la plus proche possible de la réalité”. 

Différents niveaux d’approximation et de modélisation peuvent être utilisés en fonction de la 
taille du système étudié. Plusieurs méthodes de modélisation sont utilisées dans l’optimisation 
des structures. 

Dans la dernière partie de ce chapitre, nous allons décrire les tests biologiques, ce qui donnera 


plus d'informations sur la relation structure-activité. Les résultats obtenus permettent de mettre 
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en évidence l'effet antimicrobien de la partie métallique et le groupe sulfamoyloxazolidinone 
dont l’objectif est de déterminer l’activité antimicrobienne de nos molécules étudiées. 

L’étude par docking moléculaire est une étape essentielle dans la compréhension des réactions 
biologiques entre les protéines et leurs inhibiteurs. Pour cela, le criblage virtuel vient interpréter 
les résultats obtenus et étudier théoriquement les interactions intervenant entre les composés 


synthétisés et les deux récepteurs (1IH1W et IHNJ). 
IV.2. Synthèse des ligands et des complexes 


Les ligands utilisés, dans cette étude, ont été synthétisés selon la méthode classique décrite 
dans la littérature [2]. 

Les complexes métalliques solides ont été préparés en mélangeant des solutions du ligand et 
des ions métalliques [Cu(C104)2.6H201] et [Co(C104)2.6H20] dissouts séparément dans l’éthanol 
dans un rapport ion métallique : ligand égal à 1:2. Les solutions ainsi obtenues ont été ensuite 
chauffées à reflux pendant 10 heures à une température de 50°C. Les complexes solides, formés 
à chaud, ont été filtrés et lavés plusieurs fois avec l’éthanol pour être, finalement, conservés sous 


vide en présence du déshydratant CaC12 anhydre. 


IV.3. Etude expérimentale et théorique du ligand N-(4-methoxyphenyl)-2-oxooxazolidine-3 
sulfonamide 
IV.3.1. Caractérisation de la structure du ligand N-(4-methoxyphenyl)-2-oxooxazolidine-3 


sulfonamide par DRX 


Les données du cristal pour le ligand -(4-methoxyphenyl)-2-oxooxazolidine-3 sulfonamide 
ont été obtenues avec un diffractomètre à quatre cercles Xcalibur kà géométrie automatisée 
(Oxford Diffraction, Pologne). La structure a été résolue par des méthodes directes avec SIR92 
[3] et affinée avec le raffinement de la structure des moindres carrées à matrice complète contre 
F2 à l'aide du SHELXL-2016/6 [4]. Les atomes d’hydrogène ont été ajoutés aux positions idéales 
et affinés à l’aide d’un modèle de conduite. La collecte de données, les paramètres de 
raffinement de la structure et les données cristallographiques pour le ligand sont indiqués dans le 
tableau 4.1. Le CCDC 1853764 qui contient les données cristallographiques supplémentaires 
utilisées pour les expériences peut être obtenu gratuitement du centre de données 


cristallographiques de Cambridge, à l’adresse www.ccdc.cam.ac.uk/structures 


La structure cristalline du N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide avec 


numérotation atomique est représentée sur la Figure 4.1. 
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Figure 4.1 : Structure cristalline du N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide. 
Les paramètres de raffinement de la structure et les données cristallographiques de N-(4- 


methoxyphenyl)-2-oxooxazolidine-3-sulfonamide sont rassemblés dans le Tableau 4.1. 


Tableau 4.1 : Données cristallographiquesdeN-(4-methoxyphenyl)-2-oxooxazolidine-3- 


sulfonamide 
Paramètres NÉCERORORERENDe | 
oxooxazolidine-3-sulfonamide 
Formule chimique C10 H12 N205S 
Température (K) 100(2) 
Système cristallographique Triclinique 
Groupe d’espace P21/c P-1 
a, b,c (À) 6.400(2), 8.045(2), 11.973(3) 
a,B,7(°) 98.99(2), 95.25(2), 103.32(3) 
Volume À? 587.3(3) 
Z 2 
Densité (mgm) 1.540 
F(000) 284 
0 range (°) 2.888 - 28.836 
-8h8 
Indices -10 k 10 
-15116 
Coefficient d'absorption mm! | 0.292 
Tin, Lmax 0.897, 0.988 
Diffractomètre/meas. Meth. 
No. des reflections 285 
No. des paramètres 168 
Traitement des atomes H Mixte 
R:= 0.0378 
Les indices de R ee 6 
WR2= 0.0782 


Le ligand N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide se cristallise dans le 
système de cristaux tricliniques et avec le groupe spatial P-1. Les longueurs de liaison 
sélectionnées et les angles, obtenus expérimentalement, sont donnés dans les tableaux 4.2 et 4.3. 
A titre indicatif, les distances de liaison S1-N2 et N2-H2N du groupe sulfonamide sont, 
respectivement, de 1.6032 (13) À et 0,835 (18) À. Les liaisons C1-O4 et C2-03 sont égales à 


1.2060 et 1.4593(17) À, respectivement. La valeur expérimentale de l'angle C4-N2-H2N est de 
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118.9(12) et celles des angles O1-SI-N1 (103.46(6), C7-O5-C10 (117.82(11)) et C1-03-C2 
(109.64(11)). 

L'analyse  cristallographique a démontré que La structure moléculaire de N-(4- 
methoxyphenyl)-2-oxooxazolidine-3-sulfonamideeststabilisée par des liaisons hydrogène 
intermoléculaires. À titre d'exemple, une interaction N(2)}-H(2N})e+0O(4) (2.091 À) entre le 


proton du groupement sulfonamide et l'atome d’azote O4 d'un composé voisin (Figure 4.2). 


Figure 4.2 : Vue partielle d'empilement du cristal de N-(4-methoxyphenyl)-2 oxooxazolidine-3- 


sulfonamide avec des liaisons hydrogène. 


Par modélisation moléculaire, les valeurs des longueurs de liaison et des angles au niveau 
DET dépendent, clairement, de la base et de la fonctionnelle d'échange. Sous la même base 6- 
311++G(d, p) et différents niveaux de DFT/ B3LYP, B3LYP-D3, B2PLYP et B2PLYP-D3, 
nous avons obtenu les résultats montrés dans les tableaux 4.2 et 4.3. Globalement, les valeurs 
estimées par chaque méthode ne sont pas beaucoup plus éloignées des valeurs expérimentales. 
Ainsi, la valeur expérimentale de la longueur de la liaison C4-C9 du cycle aromatique est de 
1,3930 (2) À, alors que les valeurs théoriques sont de 1,3370 et 1,4004 À en utilisant 
respectivement B3LYP et B2PLYP. Cependant, compte tenu de la dispersion de Grimme, ces 
valeurs sont considérablement améliorées. Ainsi, nous avons obtenu 1,3958 À de B3LYP-D3 et 
1,3909 À de B2PLYP-D3. De même, les valeurs expérimentales des angles C1-03-C2 et O1-S1- 
N1 égales à 109,64 (11) et 103,46 (6) sont estimées à 104,00 et 109,44 (B3LYP), 109,54 et 
105,94 (B2PLYP). Grâce aux versions de dispersion, nous avons obtenu des valeurs 
considérablement améliorées égales à 109,95 et 108,98 (B3LYP-D3), 109,62 et 102,43 
(B2PLYP-D3). 


79 


Chapitre IV Résultats et discussions 


La structure de ligand optimisée, obtenue à partir des calculs de DET, est en accord avec les 
résultats de cristallographie aux rayons X (Figure 4.1). À titre comparatif, certaines valeurs 
expérimentales et théoriques sélectionnées des longueurs de liaison ainsi que celles des angles de 
liaison de N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide sont présentées dans les 
tableaux 4.2 et 4.3. De plus, pour toutes les fonctionnelle d'échange, on peut observer que les 
valeurs optimisées sont, parfois, plus courtes ou plus longues que les celles expérimentales. Ces 
différences ont été attribuées aux effets environnementaux liés à la littérature ainsi qu’à la 
méthode DFT choisie. Par rapport aux autres groupes fonctionnels, B2PLYP-D3 permet une 
meilleure prédiction des structurelles de 


caractéristiques N-(4-methoxyphenyl)-2- 


oxooxazolidine-3-sulfonamide. 


Tableau 4.2 : Valeurs expérimentales et théoriques sélectionnées des distances de N-(4- 


methoxyphenyl)-2-oxooxazolidine-3-sulfonamide. 


Distances (À) 
re B3LYP |  +A a +A | B2PLYP | +A es +A Exp. 
S1-O1 1.4522 | 0,0284 | 1.4465 |0,0227 | 1.4609 |0,0371 | 1.4517 | 0,0279 | 1.4238(11) 
SI-N2 L.6960 0.0928 181 0,1249 1.6609 0,0577 1.6570 0,0538 | 1.6032(13) 
C1-03 1.3381 | 0,0019 | 1.3608 |0,0246 | 1.3651 |0,0289 | 1.3618 | 0,0256 | 1.3362(17) 
C2-03 1.4020 | 0,0573 | 1.4424 |0,0169 | 1.4549 |0,0044 | 1.4512 | 0,0081 | 1.4593(17) 
C7-05 1.4307 | 0,0666 | 1.3624 |0,0017 | 1.3667 | 0,0026 | 1.4219 | 0,0578 | 1.3641(17) 
CE 14680 | 0.0282 14319 0,0079 1.4338 0,006 14335 0,0063 | 1.4398(17) 
RS 10500 | 0215 10131 |0,1781 | 10151 |0,1801 | 10137 |0,1787 | 0.835(18) 
C2-H2A | 11130 | 0.123 | 10924 | 0,1024 | 1.0921 | 0.1021 | ‘10918 | 0,1018 | 0.9900 
C4-C9 13370 | 0,056 1.3958 |0,0028 | 1.4004 |0,0074 | 1.3909 | 0,0021 | 1.3930(2) 
C1-O4 1.1942 | 0,0118 | 1.1933 |0,0127 | 1.2048 |0,0012 | 1.2062 | 0,0002 |  1.2060 
CI0-HIOA | 1.1130 | 0,133 | 10952 |0,1152 } 10944 |0,1144 | 10939 ]0,1139 |  0.9800 
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Tableau 4.3 : Valeurs expérimentales et théoriques sélectionnées des angles de N-(4- 


methoxyphenyl)-2-oxooxazolidine-3-sulfonamide. 


Angles (°) 

Paramètres B3LYP +A er +A | B2PLYP | +A paca +A Exp. 
O1-S1-N2 109.50 0,08 108.71 0,87 102.24 7,34 108.54 1,04 109.58(7) 
OI-SI-N1 109.44 5,98 108.98 5,52 105.94 2,48 | 102.43 | 1,03 | 103.46(6) 

C7-O5-C10 109.50 8,32 118.21 0,39 117.89 0,07 117.79 0,03 | 117.82(11) 
C1-03-C2 104.00 5,64 109.95 0,31 109.58 0,06 | 109.62 | 0,02 | 109.64(11) 
CI1-N1-C3 111.00 0,99 110.17 1,82 109.60 2,39 109.77 2,22 | 111.99(12) 
C3-N1-SI1 124.50 1,9 118.95 7,45 118.38 8,02 118.92 7,48 | 126.40(10) 

C4-N2-H2N 120.00 1,1 116.46 2,44 115.41 3,49 114.94 3,96 118.9(12) 

N1-C3-H3B 112.29 0,49 111.33 0,47 111.20 0,6 111.37 0,43 111.80 
C9-C4-N2 120.00 1,46 119.49 1,97 119.71 1,75 120.51 0,95 | 121.46(12) 
O4-C1-03 124.27 0,62 122.89 0,76 122.08 1,57 123.20 0,45 | 123.65(11) 
O4-CI-N1 127.70 1,15 126.01 0,54 127.35 0,8 125.98 0,57 | 126.55(11) 


Pour étudier les performances des quatre fonctionnels d'échange utilisés, les coefficients de 
corrélation et la valeur de l'erreur quadratique moyenne (RMSE) entre les distances et les angles 


calculées et observées expérimentalement ont été calculés à l'aide de l'expression suivante : 


RMSE - Lier -v*) (V1) 


La représentation graphique des valeurs théoriques et expérimentales donne des coefficients 
de corrélation variant de 0.9067 à 0.9659 pour les distances et de 0.5114 à 0.8992 pour les 
angles. Il en ressort que, pour la prédiction des distances avec la B2PLYP-D3, une fonctionnelle 
double hybride (RMSE = 0.0757, R? = 0.9659) donne de meilleurs résultats et que même sans 
dispersion de B2PLYP (Fig.3), la prédiction est meilleure que celle des autres (B2PLYP, RMSE 
= 0.0794, R? = 0.9600). Il est à noter que la dispersion de Grimme améliore clairement les 
résultats du B3LYP (Tableau 4.4). 

Pour les angles, les prédictions sont moins satisfaisantes avec B3LYP (RMSE = 3.6516, R? = 
0.5114) et B2PLYP (RMSE = 3.6776, R°? = 0.7623). Néanmoins, en impliquant l'effet de la 
dispersion, les résultats théoriques deviennent en accord avec ceux obtenus expérimentalement, 


notamment avec la B2PLYP-D3 (RMSE = 2.7022, R° = 0.8992). 
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1,6 


Vocal = 0.8015xXobs +0.3026 
R2 = 0.9659 sn,” 


o 


Distances calculées A 


1,2 1,3 1,4 


Distances expérimentales (À) 


Figure 4.3 : Corrélation de distances expérimentales et calculés (B2PLYP-D3) de N- (4- 
méthoxyphényl) -2oxooxazolidine-3-sulfonamide 


Tableau 4.4 : Les coefficients de corrélation pour quatre fonctionnelles d'échange 


B3LYP B3LYP-D3 B2PLYP B2PLYP-D3 
.. | y=0.7556x +0.3647 | y = 0.8298x +0.2654 | y =0.8001x +0.3021 | y = 0.8015x +0.3026 
DistancesA R2 = 0.9067 R2 = 0.9302 R2 = 0.9644 R2 = 0.9659 
RMSE=0.0953 RMSE=0.0813 RMSE=0.0794 RMSE=0.0757 
y=0.664x+36.87 y=0.6617x+ 37.95 | y=0.7404x+ 27.39 | y = 0.7801x + 23.28 
Angles(?) R2=0.5114 R2 = 0.8312 R2 = 0.7623 R2 = 0.8992 
RMSE=3.6516 RMSE=3.0334 RMSE=3.6776 RMSE=2.7022 


IV.3.2.Etude Expérimentale et théorique des spectres infrarouge du N-(4-methoxyphenyl)- 
2-oxooxazolidine-3-sulfonamide 

La prédiction des principales fréquences de vibration FT-IR du ligand a été étudiée par 
l’utilisation de la fonctionnelle double hybride B2PLYP et la fonctionnelle hybride B3LYP. 

Les résultats des calculs, inclus dans le Tableau 4.5, sont confrontés à ceux obtenus 
expérimentalement par l'intermédiaire des représentations linéaires sur la Figure 4.4. Sur la base 
de ces constatations, examinons à titre illustratif une comparaison entre les fréquences obtenues 
B3LYP et 


théoriquement avec les 


deux fonctionnels B2PLYP et celles obtenues 
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expérimentalement pour les principaux groupements de la molécule du ligand (N-(4- 
methoxyphenyl)-2-oxooxazolidine-3-sulfonamide). 

A) Fréquences vibrationnelles du groupe N-H 

Les vibrations d’élongation de N-H donnent des bandes à 3500-3300 cm'[5] et celles de 
déformation hors plan sont indiquées à 730 cm''[6, 7]. Dans la présente étude, la bande 
d'élongation expérimentale apparait à 3283 cm ! et les valeurs prédites théoriquement sont à 3272 
cm'! et 3289 cm! par B3LYP et B2PLYP, respectivement. Les fréquences de déformation hors 
plan (yN-H) calculées sont de 732 em! (B3LYP) et de 739 cm! (B2PLYP). 

B) Fréquences vibrationnelles du groupe C-H d’un cycle aromatique 

L'existence d'un ou de plusieurs cycles aromatiques dans la structure est, normalement, facile 
à déterminer à partir des vibrations liées aux cycles C-H et C = C-C. L'élongation C-H se produit 
au-dessus de 3000 cm'! et se présente, généralement, sous la forme d'une multiplicité de bandes 
faibles à modérées par rapport à l'élongation aliphatique C-H [8]. Les bandes observées à 3120, 
3090 cm! dans le spectre expérimental FT-IR sont attribuées comme modes d'élongation C-H 
des cycles phényle. Les calculs de B3LYP donnent ces modes dans l’intervalle de 3087 à 3118 
cm'!. Le benzène possède six vibrations d’élongation, dont les quatre ayant des nombres d’ondes 
les plus élevés (se produisant vers 1600, 1580, 1490 et 1440 cm'!) sont de bonnes vibrations de 
ce groupe. En l'absence de conjugaison de cycle, la bande proche de 1580 cm! est, 
généralement, plus faible que celle de 1600 cm'!. La cinquième vibration d’élongation est active 
vers 1335 + 35 cm'!, une région qui chevauche fortement celle de la déformation dans le plan du 
C-H et dont l’intensité est, généralement, faible ou moyenne élevée [9]. Le sixième vibration 
d’élongation apparaît comme une bande faible vers 1000 cm’! dans les benzènes mono, 1,3-di et 
1,3,5-trisubstitués. Cependant, dans le benzène substitué, cette vibration de la déformation dans 
le plan est sensible aux substituants et difficile à distinguer. Pour le composé du titre, la 
déformation dans le plan du groupe phényle est observée à 1024 cm}, les B2PLYP et B3LYP 
donnent ce mode à 1020 et 1025 cm'!, respectivement. La bande observée à 1424 cm’! dans le 
spectre expérimental est attribuée aux modes vPh, avec 1420 cm! par B3LYP. Cependant, par 
B2PLYP, cette bande est observée à 1427 cm'!. 

C) Fréquences vibrationnelles des groupes C-O et C=-0 

Expérimentalement, la fréquence de vibration en mode d'élongation du motif C-O-C a été 
détectée à 1062 cm'!, valeur très proche à 1144 et 1063 cm'!, rapportées dans la littérature par 
Bhagyasree et al. [8]. En modélisation, la B2PLYP prédit avec très bon accord la fréquence de 
1064 cm'}, par contre la B3LYP donne 1059 cm. Le groupe carbonyle est contenu dans un 
grand nombre de classes différentes de composés, pour lesquels une forte bande d'absorption due 


à la vibration d’élongation C = O est observée dans l’écart de 1850-1550 cm! [7]. Pour notre 
83 


Chapitre IV Résultats et discussions 


ligand, les fréquences observées du groupe (-C = O) sont concentrées autour de 1756 cm'!. La 
B3LYP donne une valeur très proche de 1754 cm'!, alors que la B2PLYP surestime la valeur par 
9 cm. 

D) Fréquences vibrationnelles du groupe C-N 

Le classement des vibrations du C-N dans une molécule serait une entité difficile en raison de 
son occurrence en tant que mélange de plusieurs bandes dans la région. Silverstein et al. [10] ont 
attribué une absorption d’élongation C-N dans la région 1382-1266 cm! pour les amines. La 
bande de plus faible intensité à 1350 cm}, trouvée expérimentalement, a été assignée à la 
vibration symétrique (C-N). B3LYP (1347 cm!) et B2PLYP (1352 cm '!). 

E) Fréquences vibrationnelles du groupe S=0 

Normalement, les composés soufrés présentent des pics intenses dans la région 1360-1310 
cm! et 1165-1135 cem'![11] en raison des vibrations asymétriques et symétriques S = O. Pour le 
N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide, les vibrations d’élongation uS=O ont 
été observées à 1300 cm’! 


calculées sont 1306 et 1302 cm’! dans B3LYP et B2PLYP, respectivement. 


avec une très forte intensité dans le spectre FT-IR. Les valeurs 


D'une manière globale et indépendamment de la fonctionnelle utilisée, les résultats obtenus 
sont en très bon accord avec ceux obtenus expérimentalement. Ainsi, selon la Figure 4.4, le 
modèle B2PLYP semble donner les fréquences les plus précises (R2=0.9997) devant celles 
obtenues par B3LYP (R?= 0.9991). 

Le RMSE entre les fréquences vibrationnelles observées et calculées est de 23.7010 pour 
B3LYP et de 14.0432 pour B2PLYP. Cependant, la B3LYP est la fonctionnelle la plus largement 
utilisée et elle était, souvent, efficace pour prédire les fréquences de vibration des composés. 
L'utilisation du double hybride B2PLYP permet d'obtenir des résultats encore plus efficaces et 
précis dans le cas de N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide. D'autre part, le 
léger avantage, présenté par la fonctionnelle double hybride, pourrait être dû à la correction 


effectuée pour compenser les effets harmoniques. 
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Tableau 4.5 : Fréquences expérimentales et calculées (en cm!) pour certains modes de 


vibration du N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide 


Assignments B3LYP B2PLYP Exp. 
vNH 3272 3289 3283 
vC-H 3140 3124 3128 

vsCEB 2925 2923 2927 
vC=0 1754 1765 1756 
ôCH2 1530 1609 1625 
vC=C 1507 1509 1510 
vC=C 1466 1472 1488 
ôsCH3, vCH2 1380 1382 1385 
CH, vCN 1347 1352 1350 
vS=0 1306 1302 1300 
vCB, vCC 1255 1221 1234 
vCO 1059 1064 1062 
vPh, CH 1025 1020 1024 
yCE, 7 Ph 962 920 965 
yC=C 905 881 901 
yCC, , 7 Ph, yca 790 780 780 
T Ph, yC=O, yCH 758 744 750 
yC=0, yYNH 739 732 730 


*v—élongation; 5-déformation dans le plan; y-déformation hors plan; 7 — torsion; Phényle — Ph; 


as — non symétrique; s — symétrique 
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3000 y = 0,9963x + 1,3226 / 3000 y = 1,006x - 18,365 ré 
R?2 = 0,9991 Fu R2 = 0,9997 Ed 
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Figure 4.4 : Corrélation des nombres d'onde expérimentaux de N- (4méthoxyphényl) - 


2oxooxazolidine-3-sulfonamide et calculés a) (B3LYP), b) (B2PLYP) 
IV.3.3. NBOs et TD-DET de N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide 


L'analyse de la NBO a été réalisée au niveau B3LYP / 6-311++G (d, p) afin d'élucider la 
délocalisation de la densité électronique dans la molécule du titre N- (4-méthoxyphényl) - 
2oxo0oxazolidine-3-sulfonamide. 

Les valeurs E2 représentent l'énergie de stabilisation nécessaire pour diverses transitions 
possibles entre les niveaux donneur et accepteur dans la molécule. Selon les valeurs de E2, les 
transitions qui contribuent au transfert de charge intramoléculaire, par conséquent à la 
stabilisation du système, sont dues principalement au chevauchement d’orbital entre les orbitales 


liantes (x) et d'anti-liantes (x*). Ainsi, on relève : 


ñ (C14 -C15) —x * (C11-C16), 26.61 kcal / mol, 
ñ (C11-C16)—7x * (C12 -C13), 23.65 kcal / mol, 
ñ (C12-C13) —x* (C14-C15), 22.98 kcal / mol, 

Les transitions électroniques LP —x* sont impliquées entre les O1 et C2 - O10 d'une part, et 
entre O8 et C11-C16 d'autre part. Les énergies sont de 2.02 et 3.09 kcal/mol, respectivement. 

On s'attend, généralement, à ce que toutes ces transitions soient observées dans le spectre UV- 
Visible, mais certaines d'entre elles peuvent avoir une intensité très faible et peuvent ne pas être 
visibles dans le spectre. Ce qui peut également être prédit par calcul des valeurs de résistance ou 
de coefficient d'absorption de l'oscillateur pour ces trois transitions en utilisant l'analyse TD-DET 


comme suit : 
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Les transitions électroniques, au sein de la molécule dans l'éthanol comme solvant, ont été 
élucidées en utilisant la TD-DFT/B3LYP/6-311++G (d,p) level. Les énergies d'excitation 
électronique théoriques et les forces de l'oscillateur ont été calculées et les résultats ont montré 
que les bandes d'absorption les plus probables (selon les valeurs de la force d'oscillateur f) sont 
les suivantes : 228nm (f= 0.0077) and 235 nm (f=0.05). La bande apparue à 228nm correspond à 
la transition de HOMO — LUMO + 3. Le pic à 235 nm provient de la contribution de HOMO — 
LUMO + 1 ou de HOMO-2 — LUMO +3. Ces transitions ont été illustrées à la Figure 4.5. 

Ces résultats sont en très bon accord avec ceux obtenus expérimentalement par l'élargissement 


et l'extension de la bande d'absorption de 235 nm à 225 nm. 


L+1 : -0.023 eV 


, P 
à « 
%.°0 » 


H : -0.225eV H-2 : -0.294eV 


Figure 4.5 : Transitions électroniques et énergies pour N- (4-méthoxyphényl) -2- 


oxooxazolidine-3-sulfonamide 
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IV.3.4. Descripteurs de réactivité chimique globale et indices de Fukui du N- (4- 


méthoxyphényl) -2-oxooxazolidine-3-sulfonamide 


les valeurs de 7 et À de N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide sont 
calculées au niveau théorique B3LYP / 6-311++G (d,p).7 (= -Exomo) est le potentiel d'ionisation 
égale à 6.258 eV et À (= -Erumo) est l'affinité électronique de la molécule égale à 0.816 eV. 
L'énergie supérieure (HOMO) concerne la molécule la plus réactive dans les réactions avec des 
électrophiles, tandis que l'énergie inférieure (LUMO) est nécessaire pour les réactions 
moléculaires avec des nucléophiles [12]. Les valeurs calculées de GCRD telles que n, L, S, et © 
pour les molécules de N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide sont, 
respectivement, 2.721, -3.537, 0.183, 3.537 et 2.285 eV. Ces valeurs sont en faveur de la 


réactivité significative du ligand et de sa large possibilité de participer aux réactions de 


complexation en tant que donneur. Les valeurs des indices de Fukui condensées à l'atome, f 


pour l'attaque nucléophile et f, pour l'attaque électrophile, ont été calculées à l'aide de l'analyse 


de population de Hirshfeld pour les complexes étudiés. En ce qui concerne les résultats présentés 
à la Figure 4.6, on peut voir que la couleur verte se situe au voisinage des atomes de carbone 
(COS), (C14) et (C16), ce qui indique que la densité électronique augmente à ces emplacements 
sous l’attaque nucléophile. De l'autre côté, en ce qui concerne la fonction f, de Fukui, la 
couleur bleue indique que la région située à proximité des atomes d'azote (N9) et d'oxygène 


(010) perd de la densité électronique lorsque la molécule est soumise à l'attaque électrophile. 


Figure 4.6 : Représentation des indices de Fukui pour N-(4-methoxyphenyl)-2-oxooxazolidine- 


3-sulfonamide 
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IV. 4. Etude expérimental du Cu(ID et Co(IT) complexes du N-(4-methoxyphenyl)-2- 


oxooxazolidine-3-sulfonamide (L!) et N-sulfonyl imine (L?). 


IV.4.1. Formation des constantes et le diagramme de distribution des espèces des 


complexes 


Le titrage spectrophotométrique d'une solution dans l'éthanol 10°* M de Cu?* ou de Co* par 
10°M du ligand N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide (L!) est représenté sur la 
Figure 4.7. D’après les résultats spectrophotométrie, nous remarquons une apparition de deux 
nouvelles bandes, la première à 270 nm et la deuxième comme une bande large entre 310 et 330 
nm. 

Pour le titrage spectrophotométrie des mêmes ions métalliques, utilisés précédemment par N- 


sulfonyl imine (L?), nous avons observé une nouvelle bande à 340nm. 


3 


Abs 


0 I 
250 260 280 300 320 340 350 


Langueur d’onde nm 


Figure 4.7 : Le titrage spectrophotométrie d'une solution dans l'éthanol 10°* M de Cu°* 


ou de Co* par 10° M du L! 


Le traitement des données spectrophotométriques par le logiciel hypspec [13] montre 
l'existence de deux complexes ML: et ML: où M représente l'ion métallique. Les constantes de 


stability B;k sont rassemblées dans le Tableau 4.6. 
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Tableau 4.6 : Constantes de stabilité des complexes formés 


Complexes formés 


Constante de stabilité 


[Cu(LF* 


log B,, =4.65+ 0.022 


[Cu(L?»2P+ 


log, =9.27+0.010 


[Co(LF* 


log B, =3.32+ 0.033 


[Co(L')}* 


log B,, =8.02+0.064 


[Cu(L?)F* 


log f,, =2.65+0.012 


[Cu(L?»}* 


logB, =7.15+0.025 


Les valeurs des constantes de stabilité indiquent que les complexes du Cu(IT) sont plus stables 


que ceux de Co(Il) et les complexes de stockéométrie 1 :2 sont plus stables que 1 :1. Il serait 


alors instructif, à partir de ces constantes, de construire à l'aide du logiciel Hyperquad Simulation 


and Speciation (HySS) [14] les diagrammes de distribution des espèces qui montrent les 


concentrations des différentes espèces en fonction du volume du ligand . 


Le diagramme de la Figure 4.8 indique également que les ML: sont respectivement formés 


par réactions successives au dépend de ML; par déplacement des complexes. 


[Co(L'}1* 


[Co(L'}11* 


0 1.0 2.0 3.0 4.0 


Ligand volume /ml 


100 


[Cu(L'2F* 


œ 
= 


60 


40 


Formation relative to Cu?* 


20 
[Cu(L'}117#+ 


0 1.0 2.0 3.0 4.0 
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Figure 4.8 : Diagramme de distribution du titrage d'une solution de M?* (Cm?*= 10 * M) avec L! 


(CL= 10° M) 
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IV.4.2. Etude expérimentale des spectres infrarouge des complexes 


Les spectres FT-IR fournissent beaucoup d'informations précieuses sur la réaction de 
coordination et sur la formation des complexes. 

Il y a des différences significatives entre les complexes métalliques et le ligand libre. Le 
spectre FT-IR des ligands libres montre une bande nette à 3283cm !, attribuée à v (N-H), qui s'est 
déplacée à droite suite à la formation du complexe vers 3270 cm! pour les complexes de N- (4- 
méthoxyphényl) -2-oxooxazolidine-3-sulfonamide (L!) et vers 3268 cm! pour les complexes de 
N-sylfonyl imine (L?). Ceci est une caractéristique de la coordination de l'azote amino au métal. 
D'autre part, on remarque l'apparition d'une nouvelle bande dans les spectres des complexes dans 
la gamme de 488-498 cm" !, qui pourrait être attribuée à la liaison v (M-N) [15]. 

Les fréquences de vibration v (C = O) du carbonyl group dans le N- (4-méthoxyphényl) -2- 
oxooxazolidine-3-sulfonamide libre apparaissant à 1756 cm! se voient déplacer après formation 
des complexes dans l'intervalle (1760-1765cm'!). Ce qui est en faveur de la coordination du 
ligand à travers l'oxygène du groupe amide. De même, les nouvelles bandes dans la région 534- 
560 cm'!correspondent aux liaisons formées v (M-O) [16]. 

Le spectre FT-IR du complexe du cuivre à base de N-sulfonyl imine montrent une bande nette 
entre 1605 et 1613 cm'!, attribuée à v (-C=N-), qui a été diminué par rapport à une fréquence 
correspondante au ligand libre (1643 cm'!). Ceci est indicatif de la coordination de l'azote imine 
avec le centre métallique. 

Suite à la formation des complexes, les fréquences SO2 du groupement sulfonyl qui 
apparaissent aux fréquences 588 cm'!, correspondant aux modes symétriques et antisymétriques 
respectivement, restent inchangées. Ce qui indique que le soufre et les oxygènes de ce fragment 


ne participent pas à la coordination [16]. 
IV.4.3. Etude des spectres RMN du proton 


Les spectres RMN 'H des ligands N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide et 
de leurs complexes Cu (I) et Co (ID) ont également fourni des preuves supplémentaires du mode 
de liaison du ligand. Les spectres de RMN des complexes Cu (I) et Co (I) indiquaient un 
déplacement de la densité électronique du ligand au métal. Le spectre RMN 'H du complexe est 
présenté à la Figure 4.9. 

Dans les spectres de RMN 'H des complexes, les signaux du proton du groupe aminé de 
sulfonamide (-SO2-NH-) ont été observés à 10.60 ppm pour le ligand libre et à 10,84 ppm dans le 
spectre de complexe de cuivre et à 10.83 dans le spectre de complexe de cobalt en raison du NH 


non déprotoné dans le ligand. Ces différences dans les déplacements des protons, liées au 
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groupement amine dans le ligand seul et dans les complexes correspondants (A = 0.23 et 0.24 
ppm), indiquent que le mode de complexation le plus probable fait impliquer les atomes d'azote 
du ligand. 

D'autre part, dans le cycle oxazolidinone, le signal du proton les plus affecté (H2) par la 
complexation a été observé à 4.29 ppm dans le complexe du cuivre et 4.31 ppm dans celui du 
cobalt. Ce proton apparait à 4.30 ppm dans le ligand seul. Ce qui justifie l’interaction du groupe 
carbonyl de groupe oxazolidinone avec les ions métalliques. En revanche, le groupe 
méthoxyphenyl ne participe pas à la coordination, car les valeurs des déplacements des protons 


n'ont pas été affectées. 


1 
2 lo Ÿo 
H— s 
th D TS) H 384 
2.08 H 0 PA H3.9 
] Hi 1277-02 
H7 Ho #s 26 


Figure 4.9 : Spectres de RMN ‘H du complexe de cobalt 
En conclusion, tenant compte des résultats spectrophotométries du FT-IR et la RMN 'H, les 


structures des complexes sont les suivantes: 


H H . 
a H (e) H b 
CNY | _CH, le \Z ] 
D ro) D 
O & ; H : N 
L HO---CuoH  H a 
# à O . % H 
0 4 LL cc 
MC NON) De 
H H 7 NS | 
H H H 


Figure 4.10 : Les structures des complexes a :{[M(L')]* et b :[M(L?}}”*. 


IV.4.4. Analyse thermogravimétrique 


ES 


Le comportement thermique du complexe a été étudié à l'aide d'une analyse 
thermogravimétrique de la température ambiante à 800°C sous atmosphère d'azote. Sur la base 
des thermogrammes, les étapes de décomposition, les plages de température, le produit de 
décomposition ainsi que les pourcentages de perte de poids du complexe ont été examinés. La 


courbe TG et les résultats obtenus pour [Cu(L')2]?* sont présentés par la Figure 4.11. La courbe 


92 


Chapitre IV Résultats et discussions 


* se déroule en trois étapes. La 


TG montre que le processus de décomposition du [Cu(L')2]? 
première se produit entre 178 et 275 © C où le ligand éliminé entraîne une perte de masse 
expérimentale de 17,01%. Cette étape se termine par la formation de l’intermédiaire donné par la 
formule [Cu(L')]*. La deuxième est observée entre 275 et 500°C, entraînant une perte de masse 


de 48.15%. A ce stade, 2/3 moles du ligand ont été cassées et transformées en produits gazeux et 
2+ 


ES 


l'intermédiaire de formule [Cu (L')13]* s'est formé. La troisième est observée à une 
température> 500°C, entraînant une perte de masse de 29.83%. Le produit restant à 500°C 
pourrait être identifié en tant que groupe de cuivre. La réaction de décomposition thermique du 


complexe de Cu (IT) basée sur les données de TG est suggérée comme suit: 


[Cu(L ) F- 178°C-275°C [Cu(L NE 275°C-300°C [Cu(L ) T° >500 °C group de Cu 
2 1/3 
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Figure 4.11: Diagramme TG du complexe de Cu (I) (vitesse de chauffage 10°C/min, 


atmosphère N2 8.6mg). 


IV.5.Etude théorique des complexes du Cu (ID) et Co (I) a base du N-(4-methoxyphenyl)-2- 


oxooxazolidine-3-sulfonamide (L!) 


IV.5. 1. Etude de la géométrie des complexes 

La figure 4.12 montre les structures moléculaires des complexes optimisées au niveau 
B3LYP/6-311++G(d,p). Le  N-(4-methoxyphenyl)-2-oxooxazolidine-3-sulfonamide est 
coordonné à l'ion métallique (Cu ou Co(Il)) en tant que ligand bidenté. L'ion central se trouve 
dans une position où 1l est tétra coordonné par deux atomes d'azote du groupe amine et deux 
atomes d'oxygène du ligand. 

Les valeurs des longueurs de liaison et les angles sélectionnés sont regroupés dans le Tableau 
4.7. Selon les valeurs des longueurs de liaison Cu-N9 (2.090 À) et Cu-010(2.167À), l'atome de 
l'azote est donneur relativement plus fort que l'oxygène. 

Suite à la complexation, certaines longueurs de liaison du ligand sont affectées par la 
coordination avec l'ion métallique, ainsi les liaisons C20-O10 et C34-N37 sont augmentées par 
0.170 À et 0.041 À, respectivement. 

Les angles de liaison autour de l'ion métallique central (Cu) sont N15-Cu59-038 = 98.23°et 
N9- Cu59-N37 = 114.69°. Dans le ligand libre, des angles ont été calculés comme (C4- N2- 
H2N) 116.00° et (C9-C4-N2) 120.00°, ces angles sont calculés comme 109,73 et. 118.71° dans 


le complexe de cuivre, on peut dire que c'est une bonne preuve pour la formation complexe. 
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Pour Co complexe, les distances de liaison calculées autour de l'ion métallique central (Co) 
sont presque égales, Co61-N9 = 2.019À, Co61-010 = 1.941À, Co59-N37 = 2.019À et Co6l- 
038 = 1.941À. Dans le ligand libre, les distances de C1-O04 et C4-N2 sont 1.194 et 1.430 
respectivement, après la complexation ces distances sont augmentées vers 1.210 et 1.468 À, 
respectivement. 

Les angles O4-C1-03 (123.40°) et O4-CI-N1 (125.47°) sont diminués par 0.87 et 2.23° sous 


l’effet de l’interaction de l’atome d’oxygène (04) avec l’ion métallique. 


Figure 4.12 : Structures moléculaires des complexes optimisées, a)Complexe du cuivre b) 


Complexe du cobalt 
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Tableau 4.7 : Distances et angles calculés du [Cu(L!'}]*?et [Co(L')2]*? 


Distance (À) [Cu] [Co]? 
M59-N9 2.090 2.019 
M59-010 2.167 1.941 
M59-N37 2.091 2.021 
M57-038 2.167 1.941 
C2010 1.200 1.215 
O1-C2 1.330 1.323 
C11-N9 1.460 1.468 
S6- N9 1.762 1.794 
S6-07 1.447 1.447 
S6-08 1.444 1.443 
C30-N31 1.403 1.395 
C30-038 1.208 1.215 
C39-N37 1.460 1.468 
S34-035 1.447 1.447 
S34-036 1.444 1.443 
C39-C40 1.399 1.399 
C42-C43 1.400 1.401 

Angles (°) 

N9-M59-010 90.56 94.36 
010-M59-038 95.78 82.35 
N37-M59-038 90.57 94.36 
N9-M59-N37 114.43 114.69 
C12-C11-N9 118.71 119.09 
N3-C2-010 126.47 125.47 
C11-N9-M59 106.52 112.77 
S6-N9-H60 105.08 102.49 
C11-N9-H60 109,73 108.07 
029-C30-038 124.04 123.94 
010-C2-01 124.04 123.40 
010-C2-N3 126.63 125.47 
C39-N37-M59 106.51 112.77 
S34-N37-M59 116.92 117.49 


IV.5.2.Analyse de la NBO 


L’analyse de la NBO [17] et l’analyse par la théorie des perturbations du second ordre de la 


matrice de Fock fournissent des détails sur les charges naturelles, le type d’hybridation, la nature 


de la liaison et la force des interactions entre un ion métallique et les sites donneurs [18,19]. La 
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disposition électronique, la distribution électronique en 3D et la charge atomique naturelle de 


l'ion Cu (I) dans les complexes rapportés sont présentées dans le Tableau 4.8. 


Tableau 4.8 : Arrangement électronique, populations des orbitales 3d et charge de l'atome du 
métal dans les complexes étudiés. 


Complexes Arrangement électronique ds dx dy dx. d» Charge du M 


Complexe [Ar] 4S 0:33 3q C6 4p 
du Cu(D 027 5p (0.01 


1.8937 19305 19533 19042  1.9559 28.2597 


Complexe [Ar] 4$(0223d(785)4p(00) 1.9409  0.3691 1.7241 1.8786  1.9325 26.0700 
du Co(Il) 


Selon le Tableau 4.8, l'arrangement électronique des atomes de Cu et de Co dans les 
complexes étudiés est [Ar] 4S (35 3q 069 49 029 5p 00) et [Ar] 4S 022 3q (89 4p (0:01) 
respectivement. Les occupations de M 3d sont dxydx-d,zdx. dz de Cu et dxydxzdyzdx2. dzde Co. 
Pour le complexe Cu, 18 électrons sont répartis en électrons du le noyau, 10,244 électrons de 
valence (sur les atomes 4s, 3d et 4p) orbitales et 0.005 électrons en électrons de Rydberg (sur 
l’orbite 5p), ce qui donne un nombre total d’électrons de 28.259. Cela correspond à la charge 
naturelle calculée sur l'atome de Cu (+ 0.413e), qui correspond à la différence entre 28.254 et le 
nombre total d'électrons dans un atome de Cu isolé (29e). De même, le noyau, la valence, les 
électrons de Rydberg et le total des électrons sur l'atome de Co sont 18, 8.069, 0.004 et 26.070 
dans cet ordre. La charge naturelle calculée à été établie à + 0.921e. 

Le Tableau 4.9 répertorie les valeurs sélectionnées de l'énergie d'interaction de second ordre 
calculée (E2) des liaisons M-L dans les complexes étudiés. Plus les valeurs E2 sont élevées, plus 
l'interaction entre les donneurs d'électrons et les accepteurs d'électrons est intense. De plus, le 
Tableau 4.9 montre que l'énergie d'interaction entre chaque ion métallique et les atomes d'azote 
est supérieure à celle entre les ions métalliques et les atomes d'oxygène. Par contre, pour le 
complexe de cobalt, la valeur d’énergie plus élevée E2 était de 16,04 kcal/mol, ce qui correspond 
à LP (1) N9 — LP * (5). 

Les valeurs d'énergie d'interaction de second ordre calculées corroborent les résultats 
expérimentaux, tels que la stabilité des complexes de Cu (ID) et de Co (I). Les interactions 
donneur-récepteur pourraient être décrites comme un don d'électrons depuis les NBO de 
donneurs remplis d'une seule paire des atomes O et N vers les orbitales hybrides métalliques 
vides. L'analyse des orbitales naturelles des groupes de ligands indique que la deuxième paire 
isolée d'atomes de O et de N est celle qui est responsable des interactions importantes avec les 


NBO anti-liantes des ions Cu et Co. 
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Tableau 4.9 : Valeurs de l'énergie d'interaction de second ordre calculée (E2) des liaisons métal- 


L'dans les complexes étudiés 


Complexes donor— acceptor E2, kcal/mol. 

LP (1) N9 —LP*( 6) Cu 39.72 

LP (1) N37—LP*(6) Cu 49.86 
CORPS Ten LP (2) 010 — LP*(8)Cu 321 

LP (1)O 38—LP* (7) Cu 3.26 

LP (1) N9—LP*( 5) Co 16.04 

LP (1) O10—LP*( 5) Co 11.80 
COMPARE ECO LP (1) N37—LP*( 5) Co 15.85 

LP (1) 038— LP*( 5) Co 11.87 


IV.5.3. Calculs TD-DFT de complexes Cu (ID) et Co (I) 


Les longueurs d'onde d'absorption, les oscillateurs correspondants et les transitions 


électroniques verticales ont été calculés pour les complexes en milieu éthanoïque (modèle 
CPCM) par TD-DFT / B3LYP / 6-311 ++ G (d, p). Les résultats de calcul montrent que le 
complexe de cuivre a des absorptions à 274 nm (f = 0.0682), 293 nm (f = 0.0325), 312 nm (f = 


0.0621) et 331 nm (f = 0.0193). Ces résultats sont en accord avec les données expérimentales. 


Nous avons montré auparavant (Figure 4.7) que les absorptions ont été observées à 270 nm et 


dans la plage (310 - 330 nm). Les transitions verticales correspondant aux pics principaux sont: 
H — L (12%), H-2 — L + 1 (45%), H-3 — L (21%) et H — L +1 (7%). Pour le complexe de 
cobalt, les bandes d’absorption sont localisées à 271 nm (f = 0.0633), 310 nm (f = 0.0591) et 318 
nm (f = 0.0730) dont les transitions sont: H — L (10%), H -2 — L + 1 (17%), H — L + 1 (36%) 


et H— L+3(5%). 
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Tableau 4.10 : Résultats TD-DFT de complexes Cu (I) et Co (II) 


Longeur d’onde, À (nm) Force 


Expérimental TD-DFT d’oscillateurs 


Contribution 


H-2—L(17%) 
273 0.0325 H—L(38%) 
H-1—L+1 (28%) 
H—L(12%) 
270-285 281 0.0621 H-2—L+1(45%) 
H—L+1(07%) 
H-2—L+1(09%) 
292 0.0193 H-2—L+4(58%) 
H—L+4(11%) 


Complexe du 
Cu() 


H-3—L(14%) 
270-285 269 0.0591 H—L+2(30%) 
H-1—L+1 (32%) 
H—L(10%) 
H-2—L+1(17%) 
H—L+1(36%) 
H—L+3(5%) 
H-2-—L+1(6%) 
H-2-—L+4(21%) 
H—L+4(30%) 
H-1—L(2%) 


Complexe du 283 0.0730 
Co(I) 


290 0.0053 


IV.5.4. Etude théorique des spectres FT-IR des complexes 


Les fréquences vibrationnelles de la molécule sont calculées au niveau B3LYP en utilisant la 
6-311 ++ G (d, p).Les fréquences FT-IR observées pour différents modes des vibrations ont été 
présentées. La comparaison des fréquences calculées à B3LYP avec les valeurs expérimentales 
révèle la surestimation des modes vibrationnels calculés due à la négligence d'une harmonicité 
dans le système réel. L'inclusion de la corrélation d'électrons dans la théorie fonctionnelle de 
densité à certaines étendues rend les valeurs de fréquence plus petites en comparaison avec les 
données de fréquence 6-311 ++ G (d, p). Expérimentalement, les bandes à 3270 et 3272cm'! sont 
affectées aux vibrations N-H pour les complexes de cuivre et cobalt respectivement. Les valeurs 
calculées de ce mode de vibration est obtenue à 3278cm'! pour le complexe de cuivre et 3280 
cm! pour le complexe de cobalt. Pour les deux complexes, les vibrations C=O sont observées 
entre 1756-1760 cm! et les valeurs calculées de ces vibrations varient entre1769 et 1772 cm'!. 
Les vibrations S = O observées à 1310 cm! dans le spectre FT-IR ont des valeurs calculées de 


1324 cm'!. 
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IV.5.5. Potentiel électrostatique moléculaire (MEP) pour le complexe de Cu (ID) 


Le potentiel électrostatique moléculaire (MEP) pour le complexe de Cu (I) est calculé selon 
le niveau de théorie de B3LY P / 6-311++G (d, p). 

Comme le montre le MEP du complexe Cu (ID) (voir Figure 4.13), les régions ayant le 
potentiel négatif se trouvent sur les atomes électronégatifs, et les régions présentant un potentiel 
électrostatique négatif sont généralement associées à la paire isolée d'atomes électronégatifs. Le 
potentiel le plus négatif se situe sur les groupes carbonyliques, tandis que les régions les plus 
positives se situent sur les atomes d'hydrogène. L'atome de soufre entouré de couleur jaune a le 
potentiel négatif, et l'autre atome de carbone semble avoir un potentiel nul. D'après les surfaces 
MEP du complexe de Cu (Il) indiquées à la Figure 4.13, les côtés les plus réactifs sont les 


groupes carbonyle qui localisent les côtés opposés du cycle oxazolidinone. 


cs ER CE: 


Figure 4.13: Potentiel électrostatique moléculaire (MEP) du complexe de Cu (11) 
IV.5.6. Propriétés non-linéaires optiques NLO 


La première hyperpolarisabilité statique (fo) a également été calculée à l'aide de la relation 
suivante fo = (3/5) Bret est donnée dans le Tableau 4.11. 

Les valeurs calculées d'hyperpolarisabilité (Br) du complexe de Cu (ID) (3.0074x10"* esu) et 
du complexe de Co (ID) (3.4166x10"*esu) sont similaires mais trois fois plus élevées que celles 
du ligand 1.1036x10*. Habituellement, les valeurs trouvées d'hyperpolarisabilité (Bu) sont 
comparées à celles du composé de référence qui est l'urée. Cette comparaison montre que le 


ligand et en particulier les complexes ont de meilleures propriétés NLO que celles de l'urée. 
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Tableau 4.11 : Valeurs de première hyperpolarisabilité (Bt) du ligand et de leurs complexes 


Parameters EE complexe de Cu(II) Complexe de Co (II) 
Cor -102.5976 -125.1287 -15.3716 
Aa 49.0497 161.0241 325.8295 
Utot 3.4435 2.1056 3.4586 

Brotx 10*(esu) 1.1036 3.416 3.0074 


IV.6. Tests biologiques 


Un certain nombre de sulfamides, tels que le sulfanilamide, la sulfapyridine, le sulfathiazole 
et la sulfadiazine ont été les premiers agents chimio-thérapeutiques efficaces utilisés 
systématiquement pour la prévention et le traitement des infections bactériennes chez l'homme 
[20,21]. L'activité biologique prédominante de divers complexes sulfonamides métalliques a été 
rapportée avec un traitement antibactérien [22]. 

Des précédents travaux basés sur des tests biologiques, les composés oxazolidinones se sont 
montrés très efficaces pour inhiber la prolifération du parasite P. falciparumet de plusieurs 
lignées de cellules cancéreuses. Ces composé sont, donc, servi de référence pour l’évaluation des 
activités biologiques de nos analogues simplifiés [23]. 

Le potentiel antimicrobien des composés nouvellement synthétisés a été étudié vis-à-vis des 
microorganismes d'essai et exprimé comme le diamètre des zones d'inhibition selon la gélose 
méthode de diffusion sur plaque [24]. L'activité antibactérienne des composés synthétisés a été 
évalué par la méthode standard de diffusion sur disque de papier filtre. Les diamètres de zone 
d’inhibition obtenus ont été mesurés en millimètres et utilisés comme critère de l'activité 
antimicrobienne. Les solvants (DMSO) ont été inclus dans chaque expérience en tant que 
contrôles négatifs. Le DMSO a été utilisé pour dissoudre les composés testés et n'a montré 
aucune zone d'inhibition, confirmant qu'elle n'a aucune influence sur la croissance des 
microorganismes testés. 

L'activité antibactérienne du composé nouvellement synthétisé a été criblée contre certaines 
bactéries gram positives (Staphylococcus aureus (ATCC 29013)), et des bactéries gram négative 
(Escherichia Coli (ATCC27922)). Celles-ci étaient individuellement responsables de diverses 
infections et troubles, tels que l'infection de la peau, l’infection respiratoire, l’intoxication 
alimentaire, les vomissements, la diarrhée, les infections des voies urinaires, les infections des 


plaies, la fièvre typhoïde et le syndrome hémolytique et urémique [25,261]. 
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Une plaque a été préparée et le milieu stérile a été uniformément enduit sur les plaques avec 
des cultures de bactéries Gram positives et Gram négatives sélectionnées (Figure 4.14). Les 
composés étudiés imprégnés de disque de papier stérilisé placé sur une plaque de gélose sont 
ensuite inversés. Les plaques ont été incubées pendant 24 h à 37 C°. Après incubation, le résultat 
antibactérien des composés étudiés a été enregistré en tant que diamètre moyen de la zone 
d'inhibition autour du disque de papier imprégné et présenté dans le Tableau 4.12 regroupe les 


résultats de tous les tests effectués en trois répétitions. 


Gram positif Gram négatif 
Acide lipotéicroïque 
f 
Acide téichoïque Lan: Porines ?hosphclipides 
\ \ 


\\ 
Merrbrane externe | \\ 
N 


Pep:idoglycane 
_——— Espace périplasmique—"# ? 
Membrane cellulaire 7. ip. 1 
Lioprotéines 


Proteines mermbrana:res 


Figure 4.14 : Représentation schématique des parois bactériennes des 
bactéries à Gram positif et à Gram négatif 


IV.6.1. Concentration minimale inhibitrice (CMI) du composé le plus actif 


La CMI des composés synthétisés les plus puissants était déterminée par la méthode 
conventionnelle de diffusion sur disque de papier (Figure 4.15). Les bactéries ont été cultivées 
sur gélose nutritive moyenne. Les composés synthétisés les plus puissants ont été dissous dans 
du DMSO et chargés sur des disques en papier par micropipette avec différentes concentrations. 
Des disques de séchage ont été chargés à la surface de la gélose sur plaques inoculées avec 
l'organisme d'essai. L’inhibition de la croissance a été examinée après 24 h. de l'incubation à 
37°C, chaque essai a été répété trois fois. Nous avons constaté que la CMI était la plus basse 


concentration inhibant la croissance de l'organisme d'essai. 


102 


Chapitre IV Résultats et discussions 


Gradient de concentration 


( 


CMI Diamètre d’inhibition 


Figure 4.15 : Antibiogramme par diffusion en milieu gélosé 


Dans la présente étude, la sensibilité des bactéries sélectionnées vers les complexes de ligand 
et de cuivre (ID) ont été criblés par la mesure du diamètre d'inhibition. 

La valeur CMI du complexe de cuivre contre Escherichia coli a été identifiée à 3.2ug/ml. 
Une concentration élevée du complexe de cuivre était nécessaire pour l'activité anti 
Staphylococcus aureus et 50ug/ml a été considéré comme la valeur de CMI pour Staphylococcus 
aureus. Cette CMI est la petite quantité d’antibiotique (ATB) capable d’inhiber une croissance 
visible à l’œil nu. C’est l’approche la plus utilisée pour évaluer in vitro l’activité bactériostatique 
d’un antibiotique. 

Le complexe de Cu (I) a affiché une activité antibactérie nne efficace et sélective contre la 
bactérie Staphylococcus aureus. Le complexe de Cu (IT) présente activité antibactérienne plus 


importante que celle du ligand libre et les autres complexes. 


Tableau 4.12: Diamètres d’inhibition (mm) 


Composés Staphylococcus aureus Escherichia coli DMSO 
Sulfamoyoxazolidinone 2. 4 - 
Complexe de cuivre 8 8 - 
Complexe de cobalt 4 6 - 


L'activité antimicrobienne accrue des complexes peut être expliquée sur la base du concept 
Overtone [27] et la théorie de la chélation tweedy [28]. Selon le concept Overtonede, la 
perméabilité cellulaire de la membrane lipidique entourant la cellule ne favorise que le passage 
de matières liposolubles, ce qui a pour conséquence la lipo-solubilité comme étant un facteur 


important contrôlant l'activité antimicrobienne. Lors de la chélation, la polarité de l’ion 
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métallique sera réduite dans une plus grande mesure en raison du chevauchement du ligand 
partage orbital et partiel de la charge positive du métal avec les groupes de donneurs. 

Une telle complexation pourrait améliorer le caractère lipophile de l'atome de métal central, 
qui favorise par la suite sa permutation à travers les couches lipidiques de la membrane 


cellulaire, ce qui inhibe la croissance des bactéries. 
IV.7. Docking moléculaire 


Le mode d'interaction de la structure N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide 
et de ses complexes avec Escherichia Coli et Candida albicans a été présenté par les figures (16- 
19). De plus, la structure a été stabilisée par liaisons hydrogène. Les atomes donneurs et 
accepteurs importants dans les interactions de liaison hydrogène sont présentés dans le Tableau 
4.13. Dans l'interaction de liaison hydrogène, les atomes d'amine et d'oxygène agissaient comme 


un atome donneur important. 


Les résultats indiquent l’existence des interactions élevées et efficaces entre ces composés 
préparés et les deux récepteurs. La pose la plus favorable est inférieure à 2 À en RMSD (1,89 À 
pour N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide, 1,40 À pour le complexe Cu (ID 
et 1,05 pour le complexe Co (ID). Les valeurs du RMSD sont conformes aux résultats de Chikhi 
et Bensegueni (2008) [29] qui démontrent que tout programme de docking n’est performant que 


lorsque le RMSD est inférieur ou égale à 2 angstrôms. 


Les énergies de liaison les plus faibles du ligand N- (4-méthoxyphényl) -2-oxooxazolidine-3- 
sulfonamide et de ses complexes Co (I) et Cu (I) avec les récepteurs 1H1W et IHNTJ étaient de 
-6.4 et — 5.8 kcal/mol pour le ligand, -11.2 et -10.8 kcal/mol pour le complexe Co (I) et -10.2 et 
-12.4 Kkcal/mol pour le complexe Cu (I), respectivement. En outre, les résultats obtenus ont 
montré que l'énergie de liaison diminue lors de la coordination. Ainsi, les complexes Cu (D) et 
Co (I) étaient plus actifs que le ligand N- (4-méthoxyphényl) -2-oxooxazolidine-3-sulfonamide. 
Ces valeurs indiquaient que l'énergie de liaison diminue lors de la coordination. Nous pouvons 


conclure que les interactions de liaison hydrogène augmentent la bioactivité du composé. 
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Tableau 4.13 : Interactions de liaison hydrogène entre les ligands et les résidus des récepteurs 


Ligand Récepteur Nom du résidu Distance(À) 
GLUI4HN1:01 2.061 
Escherichia coli (PDB code: 1IHNJ) 
GLNISHN1:010 2.147 
Oxazolidinone ASN210:HD2103 2.132 
Candida albicans (PDB code: 1H1W) LYS111: HZIO1 1.876 
ALA162:HN O5 1.972 


GLUS57 : OE1H60 1.891 
Escherichia coli (PDB code: 1IHNJ) 
Cu() GLUS57 : OE1H61 1.929 


complex GLU166 :0E2H60 1.871 
Candida albicans (PDB code: 1H1W) 

GLU166 : OE2H61 1.742 

GLUS7 : OE1H60 1.886 

GLUS7 : OE1H61 1.958 

GLU166 : OE2H60 1.988 

GLU166 : OE2H61 1.854 


Escherichia coli (PDB code: 1IHNJ) 
Co (ID) complex 


Candida albicans (PDB code: 1H1W) 


Figure 4.16 : Structure 3D des interactions entre le ligand N-sulfamoyloxazolidinone et 1 HNJ. 


105 


Chapitre IV Résultats et discussions 


N 


fl 
fs 
S 
De 


Figure 4.18:Structure 3D des interactions entre le ligand N-sulfamoyloxazolidinone 1H1W. 


106 


Chapitre IV Résultats et discussions 


s ä 


Figure 4.19 : Structure 3D des interactions entre le complexe de cobalt et 1H1W. 
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Conclusion générale 


Le ligand ayant fait l’objet du présent travail le N- (4-méthoxyphényl) -2-oxooxazolidine-3- 


sulfonamide a été cristallisé et caractérisé expérimentalement par DRX, FT-IR et UV-visible. 


Les calculs par DET aux niveaux B3LYP, B3LYP-D3, B2LYP et B2LYP- D3 avec une 
même base 6-311++ G (d, p) ont permis de tester ces fonctionnelles quant aux prédictions des 
paramètres géométriques et des fréquences de vibration FTIR du ligand. Par ailleurs, les 
résultats sont influencés non seulement par le niveau du calcul, mais également par l’effet de 
l’incorporation de la dispersion de Grimme. Les deux fonctionnelles avec dispersion sont les 
plus efficaces dans la prédiction des longueurs de liaison, en revanche pour la prédiction des 


angles et des fréquences de vibration c’est la double hybride qui prend le dessus. 


Les analyses NBO et TD-DET ont révélé que les interactions x — x * sont l’origine de la 
stabilisation la plus dominante dans la molécule du ligand et que les transitions électroniques les 


plus probables déterminées théoriquement ont été confirmées par les résultats de UV-vis. 


Les indices de la réactivité globale et également locale déterminés ont été permis d’avoir une 
idée sur l’affinité de l’interaction de la molécule du ligand étudié avec un autre réactif et de 
prédire les sites qui seront affectés en cas d’interaction avec une autre molécule nucléophile ou 


électrophile. 


Les complexes de ce ligand avec les ions Co (I) and Cu (I) sont synthétisés et caractérisés 
par diverses techniques physico-chimiques et spectroscopiques, UV-vis, FT-IR, RMN 'H et 
ATG/DSC en solution ou à l'état solide. En solution, par spectroscopie UV-vis, l'apparition de 
bande large entre 270 et 285nm pour confirment la formation de nouvelles espèces. Le 
traitement des données spectrophotométries permet d’accéder aux valeurs des constantes de 
stabilité pour deux espèces formées ML: et ML: de stabilités très différentes. Par ailleurs les 
diagrammes de distribution des espèces pour les deux complexes ont été représentés. Les 


complexes solides ont été caractérisés par FT-IR et RMN 'H et TGA / DSC. 


La géométrie moléculaire des complexes de Cu(Il) et Co(l) ont été optimisés par les 
méthodes de la théorie fonctionnelle de la densité (DFT / B3LYP) utilisant la base 6-311 ++ G 
(d,p). Les paramètres géométriques calculés ont été comparés aux valeurs expérimentales et se 


sont révélées en très bon accord de même que pour les fréquences de vibration FT-IR. 
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Conclusion générale 
Les analyses de NBO et de TD-B3LYP montrent que l'interaction de l’ion M (I) avec 
l'atome d’azote du sulfonamide est plus forte qu’avec l’atome d’oxygène et que les longueurs 


d'onde d'excitation prédites sont en accord avec les valeurs expérimentales. 


Les propriétés NLO prédites du ligand et de ses complexes sont beaucoup plus grandes que 
celles du composé de référence. Les composés étudiés sont de bons candidats en tant que 
matériau optique non linéaire. L’hyperpolarisabilité a révélé que le complexe a les valeurs les 


plus élevées par rapport au ligand. 


L'activité antibactérienne des complexes étudiés a été évaluée in vitro, les résultats obtenus 
montrent que ces molécules sont plus actives que le ligand libre et possèdent une activité 


bactériostatique intéressante. 


La dernière partie du travail a été consacrée à l’étude théorique par le docking moléculaire de 
l’activité antibactérienne sur le ligand libre de N- (4-méthoxyphényl) -2-oxooxazolidine-3- 
sulfonamide et de son complexe Co (I) sur les récepteurs Escherichia Coli (code PDB: THNJ) et 
Candida albicans (code PDB: 1H1W). Les résultats obtenus peuvent aider à expliquer les 
interactions intermoléculaires qui influent sur l'affinité entre les inhibiteurs et la cavité de la 


protéine et de montrer que les complexes Cu (I) et Co (I) sont plus efficaces que le ligand libre. 
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